! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
! 
! The Reaction Rates File
! 
! Generated by KPP-2.2.3 symbolic chemistry Kinetics PreProcessor
!       (http://www.cs.vt.edu/~asandu/Software/KPP)
! KPP is distributed under GPL, the general public licence
!       (http://www.gnu.org/copyleft/gpl.html)
! (C) 1995-1997, V. Damian & A. Sandu, CGRER, Univ. Iowa
! (C) 1997-2005, A. Sandu, Michigan Tech, Virginia Tech
!     With important contributions from:
!        M. Damian, Villanova University, USA
!        R. Sander, Max-Planck Institute for Chemistry, Mainz, Germany
! 
! File                 : aromatics_kpp_Rates.f90
! Time                 : Fri Jun 12 02:54:05 2020
! Working directory    : /n/home08/kbates/Aromatics/MCM_update_chamber
! Equation file        : aromatics_kpp.kpp
! Output root filename : aromatics_kpp
! 
! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~



MODULE aromatics_kpp_Rates

  USE aromatics_kpp_Parameters
  USE aromatics_kpp_Global
  IMPLICIT NONE

CONTAINS



! Begin Rate Law Functions from KPP_HOME/util/UserRateLaws

!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!  User-defined Rate Law functions
!  Note: the default argument type for rate laws, as read from the equations file, is single precision
!        but all the internal calculations are performed in double precision
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

!~~~>  Arrhenius
   REAL(kind=dp) FUNCTION ARR( A0,B0,C0 )
      REAL A0,B0,C0      
      ARR =  DBLE(A0) * EXP(-DBLE(B0)/TEMP) * (TEMP/300.0_dp)**DBLE(C0)
   END FUNCTION ARR        

   REAL(kind=dp) FUNCTION TUN( A0,B0,C0 )
      REAL A0,B0,C0      
      TUN =  DBLE(A0) * EXP(DBLE(B0)/TEMP) * EXP(DBLE(C0)/TEMP**3)
   END FUNCTION TUN

   REAL(kind=dp) FUNCTION ISO1( A0,B0,C0,D0,E0,F0,G0 )
      REAL A0,B0,C0,D0,E0,F0,G0
      REAL(kind=dp) K0, K1, K2
      K0 = DBLE(D0)*EXP(DBLE(E0)/TEMP)*EXP(1.E8/TEMP**3)
      K1 = DBLE(F0)*EXP(DBLE(G0)/TEMP)
      K2 = DBLE(C0)*K0/(K0+K1)
      ISO1 =  DBLE(A0) * EXP(DBLE(B0)/TEMP) * (1.-K2)
   END FUNCTION ISO1

   REAL(kind=dp) FUNCTION ISO2( A0,B0,C0,D0,E0,F0,G0 )
      REAL A0,B0,C0,D0,E0,F0,G0
      REAL(kind=dp) K0, K1, K2
      K0 = DBLE(D0)*EXP(DBLE(E0)/TEMP)*EXP(1.E8/TEMP**3)
      K1 = DBLE(F0)*EXP(DBLE(G0)/TEMP)
      K2 = DBLE(C0)*K0/(K0+K1)
      ISO2 =  DBLE(A0) * EXP(DBLE(B0)/TEMP) * K2
   END FUNCTION ISO2

   REAL(kind=dp) FUNCTION ARRPLUS( A0,B0,C0,D0,E0 )
      REAL A0,B0,C0,D0,E0 
      REAL(kind=dp) K0     
      K0 = DBLE(A0) * (DBLE(D0)+(TEMP*DBLE(E0)))
      ARRPLUS =  K0 * EXP(-DBLE(B0)/TEMP) * (TEMP/300.0_dp)**DBLE(C0)
      IF (ARRPLUS < 0.0) ARRPLUS = 0.0
   END FUNCTION ARRPLUS

   REAL(kind=dp) FUNCTION TUNPLUS( A0,B0,C0,D0,E0 )
      REAL A0,B0,C0,D0,E0
      REAL(kind=dp) K0
      K0 = DBLE(A0) * (DBLE(D0)+(TEMP*DBLE(E0)))
      TUNPLUS = K0 * EXP(DBLE(B0)/TEMP) * EXP(DBLE(C0)/TEMP**3)
      IF (TUNPLUS < 0.0) TUNPLUS = 0.0
   END FUNCTION TUNPLUS

!~~~> Simplified Arrhenius, with two arguments
!~~~> Note: The argument B0 has a changed sign when compared to ARR
   REAL(kind=dp) FUNCTION ARR2( A0,B0 )
      REAL A0,B0           
      ARR2 =  DBLE(A0) * EXP( DBLE(B0)/TEMP )              
   END FUNCTION ARR2          

   REAL(kind=dp) FUNCTION EP2(A0,C0,A2,C2,A3,C3)
      REAL A0,C0,A2,C2,A3,C3
      REAL(kind=dp) K0,K2,K3            
      K0 = DBLE(A0) * EXP(-DBLE(C0)/TEMP)
      K2 = DBLE(A2) * EXP(-DBLE(C2)/TEMP)
      K3 = DBLE(A3) * EXP(-DBLE(C3)/TEMP)
      K3 = K3*CFACTOR*1.0E6_dp
      EP2 = K0 + K3/(1.0_dp+K3/K2 )
   END FUNCTION EP2

   REAL(kind=dp) FUNCTION EP3(A1,C1,A2,C2) 
      REAL A1, C1, A2, C2
      REAL(kind=dp) K1, K2      
      K1 = DBLE(A1) * EXP(-DBLE(C1)/TEMP)
      K2 = DBLE(A2) * EXP(-DBLE(C2)/TEMP)
      EP3 = K1 + K2*(1.0E6_dp*CFACTOR)
   END FUNCTION EP3 

   REAL(kind=dp) FUNCTION EPO(A1,E1,M1)
      REAL A1, E1, M1
      REAL(kind=dp) K1      
      K1 = 1.0_dp/(DBLE(M1) * CFACTOR + 1.0_dp)
      EPO = DBLE(A1) * EXP(DBLE(E1)/TEMP) *  K1
   END FUNCTION EPO

   REAL(kind=dp) FUNCTION KCO(A1,M1)
      REAL A1, M1
      KCO = DBLE(A1) * (1.0_dp + (CFACTOR / DBLE(M1)))
   END FUNCTION KCO

   REAL(kind=dp) FUNCTION FALL ( A0,B0,C0,A1,B1,C1,CF)
      REAL A0,B0,C0,A1,B1,C1,CF
      REAL(kind=dp) K0, K1     
      K0 = DBLE(A0) * EXP(-DBLE(B0)/TEMP)* (TEMP/300.0_dp)**DBLE(C0)
      K1 = DBLE(A1) * EXP(-DBLE(B1)/TEMP)* (TEMP/300.0_dp)**DBLE(C1)
      K0 = K0*CFACTOR*1.0E6_dp
      K1 = K0/K1
      FALL = (K0/(1.0_dp+K1))*   &
           DBLE(CF)**(1.0_dp/(1.0_dp+(LOG10(K1))**2))
   END FUNCTION FALL

   REAL(kind=dp) FUNCTION PAN ( A0,B0,C0,A1,B1,C1,CF)
      REAL A0,B0,C0,A1,B1,C1,CF
      REAL(kind=dp) K0, K1, KR, NC, F     
      K0 = DBLE(A0) * EXP(DBLE(B0)/TEMP) * (TEMP/300.0_dp)**DBLE(C0)
      K1 = DBLE(A1) * EXP(DBLE(B1)/TEMP) * (TEMP/300.0_dp)**DBLE(C1)
      K0 = K0*CFACTOR
      KR = K0/K1
      NC = 0.75_dp-1.27_dp*(LOG10(DBLE(CF)))
      F  = 10.0_dp**(LOG10(DBLE(CF))/(1+(LOG10(KR)/NC)**2))
      PAN = K0*K1*F/(K0+K1)
   END FUNCTION PAN

   REAL(kind=dp) FUNCTION VILA ( A0,B0,C0,A1,B1,C1)
      REAL A0,B0,C0,A1,B1,C1
      REAL(kind=dp) K0, K1     
      K0 = DBLE(A0) * EXP(DBLE(B0)/TEMP) * (TEMP/300.0_dp)**DBLE(C0)
      K1 = DBLE(A1) * EXP(DBLE(B1)/TEMP) * (TEMP/300.0_dp)**DBLE(C1)
      K0 = K0*CFACTOR
      VILA = K0*K1*0.35_dp/(K0+K1)
   END FUNCTION VILA

   REAL(kind=dp) FUNCTION NIT ( A0,B0,C0,n,X0,Y0)
      REAL A0,B0,C0,n,X0,Y0
      REAL(kind=dp) K0, K1, K2, K3, K4     
      K0 = 2.0E-22_dp * EXP(DBLE(n))
      K1 = 4.3E-1_dp*(TEMP/298.0_dp)**(-8)
      K0 = K0*CFACTOR
      K1 = K0/K1
      K2 = (K0/(1.0_dp+K1))*   &
           4.1E-1_dp**(1.0_dp/(1.0_dp+(LOG10(K1))**2))
      K3 = K2/(K2+DBLE(C0))
      K4 = DBLE(A0)*(DBLE(X0)-TEMP*DBLE(Y0))
      NIT = K4 * EXP(DBLE(B0)/TEMP) * K3
   END FUNCTION NIT

   REAL(kind=dp) FUNCTION ALK ( A0,B0,C0,n,X0,Y0)
      REAL A0,B0,C0,n,X0,Y0
      REAL(kind=dp) K0, K1, K2, K3, K4  
      K0 = 2.0E-22_dp * EXP(DBLE(n))
      K1 = 4.3E-1_dp*(TEMP/298.0_dp)**(-8)
      K0 = K0*CFACTOR
      K1 = K0/K1
      K2 = (K0/(1.0_dp+K1))*   &
           4.1E-1_dp**(1.0_dp/(1.0_dp+(LOG10(K1))**2))
      K3 = DBLE(C0)/(K2+DBLE(C0))
      K4 = DBLE(A0)*(DBLE(X0)-TEMP*DBLE(Y0))
      ALK = K4 * EXP(DBLE(B0)/TEMP) * K3
   END FUNCTION ALK

   REAL(kind=dp) FUNCTION EXTARR ( A0,B0,C0 )
      REAL A0,B0,C0      
      EXTARR =  DBLE(A0) * EXP(DBLE(C0)/TEMP) * (TEMP/300.0_dp)**DBLE(B0)
   END FUNCTION EXTARR

   REAL(kind=dp) FUNCTION IEMISS ( A0,B0,C0 )
      REAL A0,B0,C0 
      REAL(kind=dp) CT, CD, CL, BS, AT
      AT = DBLE(A0) * 8.314_dp * TEMP
      BS = DBLE(B0) * SUN
      CL = 2.88E-3_dp*BS/((1.0_dp+7.29E-6_dp*BS**2.0_dp)**0.5_dp)
      CT = EXP(9.5E+4_dp*(TEMP-DBLE(A0))/AT)
      CD = 1 + EXP(2.3E+5_dp*(TEMP-3.14E+2_dp)/AT)
      IEMISS = DBLE(C0) * CT * CL / CD
   END FUNCTION IEMISS

   REAL(kind=dp) FUNCTION MCMHV ( A0,B0,C0 )
      REAL A0,B0,C0
      IF (SUN.GT.0) THEN
         MCMHV = DBLE(A0) * SUN**DBLE(B0) * EXP(-DBLE(C0)/SUN)
      ELSE
         MCMHV = 0.0
      ENDIF
   END FUNCTION MCMHV

  !---------------------------------------------------------------------------

  ELEMENTAL REAL(kind=dp) FUNCTION k_3rd(temp,cair,k0_300K,n,kinf_300K,m,fc)

    INTRINSIC LOG10

    REAL(kind=dp), INTENT(IN) :: temp      ! temperature [K]
    REAL(kind=dp), INTENT(IN) :: cair      ! air concentration [molecules/cm3]
    REAL, INTENT(IN) :: k0_300K   ! low pressure limit at 300 K
    REAL, INTENT(IN) :: n         ! exponent for low pressure limit
    REAL, INTENT(IN) :: kinf_300K ! high pressure limit at 300 K
    REAL, INTENT(IN) :: m         ! exponent for high pressure limit
    REAL, INTENT(IN) :: fc        ! broadening factor (usually fc=0.6)
    REAL(kind=dp) :: zt_help, k0_T, kinf_T, k_ratio

    zt_help = 300._dp/temp
    k0_T    = k0_300K   * zt_help**(n) * cair ! k_0   at current T
    kinf_T  = kinf_300K * zt_help**(m)        ! k_inf at current T
    k_ratio = k0_T/kinf_T
    k_3rd   = k0_T/(1._dp+k_ratio)*fc**(1._dp/(1._dp+LOG10(k_ratio)**2))

  END FUNCTION k_3rd

  !---------------------------------------------------------------------------

  ELEMENTAL REAL(kind=dp) FUNCTION k_arr (k_298,tdep,temp)
    ! Arrhenius function

    REAL,     INTENT(IN) :: k_298 ! k at T = 298.15K
    REAL,     INTENT(IN) :: tdep  ! temperature dependence
    REAL(kind=dp), INTENT(IN) :: temp  ! temperature

    INTRINSIC EXP

    k_arr = k_298 * EXP(tdep*(1._dp/temp-3.3540E-3_dp)) ! 1/298.15=3.3540e-3

  END FUNCTION k_arr

  REAL(kind=dp) FUNCTION OH_O1D (J, H2O, TEMP, CFACTOR)
  REAL*8 J, H2O, TEMP, CFACTOR
  REAL*8 K1, K2, K3
  REAL*8 N2, O2

  N2=0.79d0
  O2=0.21d0

  K1=1.63d-10*exp(60d0/TEMP)
  K2=2.15d-11*exp(110d0/TEMP)
  K3=3.30d-11*exp(55d0/TEMP)

  OH_O1D=J*K1*H2O/(K1*H2O+K2*N2*CFACTOR+K3*O2*CFACTOR)
  END FUNCTION OH_O1D

  REAL(kind=dp) FUNCTION HO2_H2O(H2O,TEMP)
  REAL*8 TEMP, H2O
  HO2_H2O=1+1.4D-21*H2O*EXP(2200/TEMP)
  END FUNCTION HO2_H2O

  REAL(kind=dp) FUNCTION GCARR( A0,B0,C0 )
      REAL A0,B0,C0
      GCARR =  DBLE(A0) * EXP(DBLE(C0)/TEMP) * (300._dp/TEMP)**DBLE(B0)
  END FUNCTION GCARR

  REAL(kind=dp) FUNCTION GC_HO2NO3( A0,B0,C0,A1,B1,C1 )
      REAL A0,B0,C0,A1,B1,C1
      REAL(kind=dp) :: R0,R1
      R0 =  DBLE(A0) * EXP(DBLE(C0)/TEMP) * (300._dp/TEMP)**DBLE(B0)
      R1 =  DBLE(A1) * EXP(DBLE(C1)/TEMP) * (300._dp/TEMP)**DBLE(B1)

      GC_HO2NO3 = (R0+R1*CFACTOR)*(1.D0+1.4E-21_dp*CFACTOR*1.0E-2_dp &
                   * EXP(2200.E+0_dp/TEMP))
  END FUNCTION GC_HO2NO3

  REAL(kind=dp) FUNCTION GC_TBRANCH( A0,B0,C0,A1,B1,C1 )
! Temperature Dependent Branching Ratio                                                                                                     
      REAL A0,B0,C0,A1,B1,C1
      REAL(kind=dp) :: R0,R1
      R0 =  DBLE(A0) * EXP(DBLE(C0)/TEMP) * (300._dp/TEMP)**DBLE(B0)
      R1 =  DBLE(A1) * EXP(DBLE(C1)/TEMP) * (300._dp/TEMP)**DBLE(B1)

      GC_TBRANCH = R0/(1.D0+R1)

  END FUNCTION GC_TBRANCH

  REAL(kind=dp) FUNCTION GC_RO2HO2( A0,B0,C0,A1,B1,C1 )
! Carbon Dependence of RO2+HO2                                                                                                              
      REAL A0,B0,C0,A1,B1,C1
      REAL(kind=dp) :: R0,R1
      R0 =  DBLE(A0) * EXP(DBLE(C0)/TEMP) * (300._dp/TEMP)**DBLE(B0)
      R1 =  DBLE(A1) * EXP(DBLE(C1)/TEMP) * (300._dp/TEMP)**DBLE(B1)

      GC_RO2HO2 = R0*(1E0-EXP(-0.245E0*R1))

  END FUNCTION GC_RO2HO2

  REAL(kind=dp) FUNCTION GC_DMSOH( A0,B0,C0,A1,B1,C1 )
! DMS+OH+O2                                                                                                                                 
    REAL A0,B0,C0,A1,B1,C1
    REAL(kind=dp) :: R0,R1
    R0 =  DBLE(A0) * EXP(DBLE(C0)/TEMP) * (300._dp/TEMP)**DBLE(B0)
    R1 =  DBLE(A1) * EXP(DBLE(C1)/TEMP) * (300._dp/TEMP)**DBLE(B1)
!    GC_DMSOH = R0/(1e0_dp+R1*0.2095e0_dp)                                                                                                  
    GC_DMSOH = (R0*CFACTOR*0.2095e0_dp)/(1e0_dp+R1*0.2095e0_dp)

  END FUNCTION GC_DMSOH

  REAL(kind=dp) FUNCTION GC_GLYXNO3( A0,B0,C0 )
! ---  K = K1*([O2]+3.5D18)/(2*[O2]+3.5D18) --- HO2+2*CO branch of GLYX+OH/NO3                                                              
    REAL A0,B0,C0
    REAL(kind=dp) R0
    REAL(kind=dp) O2

    O2 = CFACTOR*0.2095e0_dp
    R0 =  DBLE(A0) * EXP(DBLE(C0)/TEMP) * (300._dp/TEMP)**DBLE(B0)
    GC_GLYXNO3 = R0*(O2+3.5E+18_dp)/(2.E+0_dp*O2+3.5E+18_dp)

  END FUNCTION GC_GLYXNO3

  REAL(kind=dp) FUNCTION GC_OHHNO3( A0,B0,C0,A1,B1,C1,A2,B2,C2 )
! ---  OH + HNO3:   K = K0 + K3[M] / (1 + K3[M]/K2)  ------                                                                                 
    REAL A0,B0,C0,A1,B1,C1,A2,B2,C2
    REAL(kind=dp) R0,R1,R2
    R0 =  DBLE(A0) * EXP(DBLE(C0)/TEMP) * (300._dp/TEMP)**DBLE(B0)
    R1 =  DBLE(A1) * EXP(DBLE(C1)/TEMP) * (300._dp/TEMP)**DBLE(B1)
    R2 =  CFACTOR*(DBLE(A2) * EXP(DBLE(C2)/TEMP) * (300._dp/TEMP)**DBLE(B2))
    GC_OHHNO3 = R0 + R2/(1.E0_dp + R2/R1)

  END FUNCTION GC_OHHNO3

  REAL(kind=dp) FUNCTION GC_GLYCOHA( A0,B0,C0 )
!                                                                                                                                           
    REAL A0,B0,C0,R0,GLYC_FRAC
    R0 =  DBLE(A0) * EXP(DBLE(C0)/TEMP) * (300._dp/TEMP)**DBLE(B0)
    GLYC_FRAC=1e+0_dp-11.0729e+0_dp*EXP(-(1._dp/73._dp)*TEMP)
    IF (GLYC_FRAC<0e+0_dp) GLYC_FRAC=0e+0_dp
    GC_GLYCOHA=R0*GLYC_FRAC

  END FUNCTION GC_GLYCOHA

  REAL(kind=dp) FUNCTION GC_GLYCOHB( A0,B0,C0 )
!                                                                                                                                           
    REAL A0,B0,C0
    REAL(kind=dp) :: R0,GLYC_FRAC
    R0 =  DBLE(A0) * EXP(DBLE(C0)/TEMP) * (300._dp/TEMP)**DBLE(B0)
    GLYC_FRAC=1e+0_dp-11.0729e+0_dp*EXP(-(1._dp/73._dp)*TEMP)
    IF (GLYC_FRAC<0e+0_dp) GLYC_FRAC=0e+0_dp
    GC_GLYCOHB=R0*(1e0_dp-GLYC_FRAC)

  END FUNCTION GC_GLYCOHB

  REAL(kind=dp) FUNCTION GC_HACOHA( A0,B0,C0 )
!                                                                                                                                           
    REAL A0,B0,C0
    REAL(kind=dp) :: R0,HAC_FRAC
    R0 =  DBLE(A0) * EXP(DBLE(C0)/TEMP) * (300._dp/TEMP)**DBLE(B0)
    HAC_FRAC=1e+0_dp-23.7e+0_dp*EXP(-(1._dp/60._dp)*TEMP)
    IF (HAC_FRAC<0e+0_dp) HAC_FRAC=0e+0_dp
    GC_HACOHA=R0*HAC_FRAC

  END FUNCTION GC_HACOHA

  REAL(kind=dp) FUNCTION GC_HACOHB( A0,B0,C0 )
!                                                                                                                                           
    REAL A0,B0,C0
    REAL(kind=dp) :: R0,HAC_FRAC
    R0 =  DBLE(A0) * EXP(DBLE(C0)/TEMP) * (300._dp/TEMP)**DBLE(B0)
    HAC_FRAC=1e+0_dp-23.7e+0_dp*EXP(-(1._dp/60._dp)*TEMP)
    IF (HAC_FRAC<0e+0_dp) HAC_FRAC=0e+0_dp
    GC_HACOHB=R0*(1.E0_dp-HAC_FRAC)

  END FUNCTION GC_HACOHB

  REAL(kind=dp) FUNCTION GC_OHCO( A0,B0,C0 )

    REAL A0,B0,C0,R0
    REAL KLO1,KLO2,KHI1,KHI2,XYRAT1,XYRAT2,BLOG1,BLOG2,FEXP1,FEXP2
    REAL KCO1,KCO2,KCO

    R0 =  DBLE(A0) * EXP(DBLE(C0)/TEMP) * (300._dp/TEMP)**DBLE(B0)
    R0 = R0 * (1.E+0_dp + 0.6e+0_dp*9.871E7_dp)

    ! new OH+CO rate from JPL2006.                                                                                                          
    KLO1=5.9E-33_dp*(300._dp/TEMP)**(1.E+0_dp)
    KHI1=1.1E-12_dp*(300._dp/TEMP)**(-1.3E0_dp)
    XYRAT1=KLO1*CFACTOR/KHI1
    BLOG1=LOG10(XYRAT1)
    FEXP1=1.E+0_dp/(1.E+0_dp+BLOG1*BLOG1)
    KCO1=KLO1*CFACTOR*0.6**FEXP1/(1.e+0_dp+XYRAT1)
    KLO2=1.5E-13_dp*(300._dp/TEMP)**(0.E+0_dp)
    KHI2=2.1e+09_dp *(300._dp/TEMP)**(-6.1E+0_dp)
    XYRAT2=KLO2*CFACTOR/KHI2
    BLOG2=LOG10(XYRAT2)
    FEXP2=1.E+0_dp/(1.E+0_dp+BLOG2*BLOG2)
    KCO2=KLO2*0.6**FEXP2/(1.e+0_dp+XYRAT2)
    KCO=KCO1+KCO2
    GC_OHCO=KCO

  END FUNCTION GC_OHCO

  REAL(kind=dp) FUNCTION GC_RO2NO( B,A0,B0,C0,A1,B1,C1 )
    ! ---  K = K1*(1-FYRNO3(K2,M,T))  ---  abstraction branch of RO2+NO                                                                     
    CHARACTER(*) B !Branch Toggle                                                                                                           
    REAL A0,B0,C0,A1,B1,C1
    REAL(kind=dp) :: R0,R1
    REAL(kind=dp) :: YYYN, XXYN,  AAA,  RARB, ZZYN
    REAL(kind=dp) :: XF, ALPHA, Y300, BETA, XMINF, XM0
    REAL(kind=dp) :: FYRNO3
    Y300 = 0.826
    ALPHA= 1.94e-22
    BETA = 0.97
    XM0  = 0.
    XMINF= 8.1
    XF   = 0.411

    R0 =  DBLE(A0) * EXP(DBLE(C0)/TEMP) * (300._dp/TEMP)**DBLE(B0)
    R1 =  DBLE(A1) * EXP(DBLE(C1)/TEMP) * (300._dp/TEMP)**DBLE(B1)

    ! Initialize static variables                                                                                                           

    XXYN   = ALPHA*EXP(BETA*R1)*CFACTOR*((300./TEMP)**XM0)
    YYYN   = Y300*((300./TEMP)**XMINF)
    AAA    = LOG10(XXYN/YYYN)
    ZZYN   = 1./(1.+ AAA*AAA )
    RARB   = (XXYN/(1.+ (XXYN/YYYN)))*(XF**ZZYN)
    FYRNO3 = RARB/(1. + RARB)
    IF (trim(B) .eq. 'A') THEN
       GC_RO2NO     = R0 * FYRNO3
    ELSEIF (trim(B) .eq. 'B') THEN
       GC_RO2NO     = R0 * (1.E+0_dp - FYRNO3)
    ENDIF

  END FUNCTION GC_RO2NO

  REAL(kind=dp) FUNCTION GCJPL3( k0_300, n, ki_300, m)
!  Functions given in JPL Booklet                                                                                                           
      REAL k0_300, n, ki_300,m
      REAL k0, ki

      k0=k0_300*((TEMP/300.d0)**(-n))
      ki=ki_300*((TEMP/300.d0)**(-m))

!      GCJPL3=(k0*CFACTOR)/(1+k0*CFACTOR/ki)*0.6** &                                                                                          
!       ((1+((LOG10(k0*CFACTOR/ki))**2d0)**-1.0d0))                                                                                          
      GCJPL3=(k0/(1.d0+k0/(ki/CFACTOR)))*0.6** &
        ((1+((log10(k0/(ki/CFACTOR)))**2d0)**1.0e0))
      GCJPL3=GCJPL3*CFACTOR
  END FUNCTION GCJPL3

  REAL(kind=dp) FUNCTION GCJPLEQ( A0,B0,C0,A1,B1,C1,A2,B2,C2,FV,FCT1,FCT2 )
 ! Function calculates the rate constant of the forward reaction                                                                            
 ! calculates the equilibrium constant                                                                                                      
 ! Find the backwards reaction by K=kforward/kbackwards                                                                                     
       REAL A0,B0,C0,A1,B1,C1
       REAL(kind=dp) :: R0,R1
       REAL, OPTIONAL :: A2,B2,C2,FV,FCT1,FCT2 !If a P-dependent rxn                                                                        

       ! Calculate Backwards reaction                                                                                                       
       R0 = GCARR( A0,B0,C0 )

       ! Calculate forwards reaction                                                                                                        
       IF (present(A2)) THEN ! P-dependent                                                                                                  
          IF (present(B2) .and. present(C2) .and. present(FV) &
                          .and. present(FCT1) .and. present(FCT2)) THEN
             R1 = GCJPLPR( A1,B1,C1,A2,B2,C2,FV,FCT1,FCT2)
          ELSE ! Missing params!                                                                                                            
             write(*,'(a)') 'GCJPLEQ: Missing parameters for P-dependent reaction.'
             write(*,'(a)') 'GCJPLEQ: Returning zero'
             GCJPLEQ = 0.E0
             RETURN
          ENDIF
       ELSE
          R1 = gcarr( A1,B1,C1 ) !Std. Arrhenius eqn.                                                                                       
       ENDIF

       GCJPLEQ=R1/R0
  END FUNCTION GCJPLEQ

  REAL(kind=dp) FUNCTION GCJPLPR(A0,B0,C0,A1,B1,C1,FV,FCT1,FCT2)
! * PRESSURE-DEPENDENT EFFECTS                                                                                                              
! * ADD THE THIRD BODY EFFECT FOR PRESSURE DEPENDENCE OF RATE                                                                               
! * COEFFICIENTS.                                                                                                                           
! A0 B0, & C0 are the Arrhenius parameters for the lower-limit                                                                              
! rate. A1, B1 & C1 are the upper-limit parameters.                                                                                         
! FV is the falloff curve paramter, (SEE ATKINSON ET. AL (1992)                                                                             
! J. PHYS. CHEM. REF. DATA 21, P. 1145). USUALLY = 0.6                                                                                      
!                                                                                                                                           
       REAL A0,B0,C0,A1,B1,C1,FV,FCT1,FCT2
       REAL FCT,XYRAT,BLOG,RLOW,RHIGH,FEXP

       RLOW  = GCARR( A0,B0,C0 )*CFACTOR
       RHIGH = GCARR( A1,B1,C1 )

       IF     (FCT2.NE.0.) THEN
             FCT            = EXP(-TEMP / FCT1) + EXP(-FCT2 / TEMP)
             XYRAT          = RLOW/RHIGH
             BLOG           = LOG10(XYRAT)
             FEXP           = 1.e+0_dp / (1.e+0_dp + BLOG * BLOG)
             GCJPLPR        = RLOW*FCT**FEXP/(1e+0_dp+XYRAT)
       ELSEIF (FCT1.NE.0.) THEN
             FCT            = EXP(-TEMP / FCT1)
             XYRAT          = RLOW/RHIGH
             BLOG           = LOG10(XYRAT)
             FEXP           = 1.e+0_dp / (1.e+0_dp + BLOG * BLOG)
             GCJPLPR        = RLOW*FCT**FEXP/(1e+0_dp+XYRAT)
       ELSE
             XYRAT          = RLOW/RHIGH
             BLOG           = LOG10(XYRAT)
             FEXP           = 1.e+0_dp / (1.e+0_dp + BLOG * BLOG)
             GCJPLPR        = RLOW*FV**FEXP/(1e+0_dp+XYRAT)
       ENDIF


  END FUNCTION GCJPLPR

  REAL(kind=dp) FUNCTION GCIUPAC3(ko_300,n,ki_300,m,Fc)
! Function calcualtes the rate constant of 3 body reaction using IUPAC                                                                      
! methology                                                                                                                                 
  REAL ko_300,n,ki_300,m,Fc
  REAL ko, ki, F, NN

  ko=ko_300*((TEMP/300.e0)**n)*CFACTOR
  ki=ki_300*((TEMP/300.e0)**m)

  NN=0.75-1.27*LOG10(Fc)
  F=10.0**(LOG10(Fc)/(1.0e0+(LOG10(ko/ki)/NN)**2.0))


  GCIUPAC3=ko/(1+ko/ki)*F
  END FUNCTION GCIUPAC3


!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!  End of User-defined Rate Law functions
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

! End Rate Law Functions from KPP_HOME/util/UserRateLaws


! Begin INLINED Rate Law Functions


! End INLINED Rate Law Functions

! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
! 
! Update_SUN - update SUN light using TIME
!   Arguments :
! 
! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

  SUBROUTINE Update_SUN()
    SUN = 5.0e+1_dp/1.0e+1_dp
 END SUBROUTINE Update_SUN

! End of Update_SUN function
! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~


! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
! 
! Update_RCONST - function to update rate constants
!   Arguments :
! 
! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

SUBROUTINE Update_RCONST ( )




! Begin INLINED RCONST

 

 RO2 = & 
 C(ind_OXYLO2) + C(ind_OXYBIPERO2) + C(ind_MXYLO2) + C(ind_MXYBIPERO2) + & 
   C(ind_PXYLO2) + C(ind_PXYBIPERO2) + C(ind_BZBIPERO2) + &  
   C(ind_C6H5CH2O2) + C(ind_TLBIPERO2) + C(ind_OXYMUCCO3) + &  
   C(ind_OXYMUCO2) + C(ind_MC6CO2O2) + C(ind_NOXYOLO2) + C(ind_OXYOLO2) + &  
   C(ind_MXYMUCCO3) + C(ind_MXYMUCO2) + C(ind_C726CO5O2) + C(ind_CH3CO3) &  
   + C(ind_MXYOLO2) + C(ind_NMXYOLO2) + C(ind_PXYMUCCO3) + &  
   C(ind_PXYMUCO2) + C(ind_C6M5CO2O2) + C(ind_NPXYOLO2) + C(ind_PXYOLO2) &  
   + C(ind_BZEMUCCO3) + C(ind_BZEMUCO2) + C(ind_C5DIALO2) + &  
   C(ind_NPHENO2) + C(ind_PHENO2) + C(ind_CRESO2) + C(ind_NCRESO2) + &  
   C(ind_TLEMUCCO3) + C(ind_TLEMUCO2) + C(ind_C615CO2O2) + C(ind_OXYLCO3) &  
   + C(ind_OXYL1O2) + C(ind_CH3O2) + C(ind_EPXM2DLCO3) + C(ind_C4MCO2O2) &  
   + C(ind_DM123O2) + C(ind_MXYLCO3) + C(ind_MXYL1O2) + C(ind_C3MCODBCO3) &  
   + C(ind_C5CO14O2) + C(ind_EPXMDLCO3) + C(ind_C3MDIALO2) + &  
   C(ind_MXY1O2) + C(ind_PXYLCO3) + C(ind_PXYL1O2) + C(ind_PXY1O2) + &  
   C(ind_MALDIALCO3) + C(ind_EPXDLCO3) + C(ind_C3DIALO2) + &  
   C(ind_MALDIALO2) + C(ind_C6H5O2) + C(ind_C6H5CO3) + C(ind_NBZFUO2) + &  
   C(ind_BZFUO2) + C(ind_BIACETO2) + C(ind_C5DICARBO2) + &  
   C(ind_MC4CODBCO3) + C(ind_MC5CO2OHO2) + C(ind_NOXYFUO2) + &  
   C(ind_C6OTKETO2) + C(ind_OXYFUO2) + C(ind_NTLFUO2) + C(ind_TLFUO2) + &  
   C(ind_OXCATEC1O2) + C(ind_C4MCODBCO3) + C(ind_C5MCO2OHO2) + &  
   C(ind_NMXYFUO2) + C(ind_C23O3MO2) + C(ind_MXYFUO2) + C(ind_MC3CODBCO3) &  
   + C(ind_C4M2ALOHO2) + C(ind_NPXYFUO2) + C(ind_MCOCOMOXO2) + &  
   C(ind_PXYFUO2) + C(ind_HCOCH2O2) + C(ind_MXCATEC1O2) + C(ind_DMKOHO2) &  
   + C(ind_C4CO2O2) + C(ind_PXCATEC1O2) + C(ind_HCOCOHCO3) + &  
   C(ind_CATEC1O2) + C(ind_MCATEC1O2) + C(ind_NTMB1FUO2) + &  
   C(ind_TMB1FUO2) + C(ind_HCOCO3) + C(ind_MECOACETO2) + C(ind_NOXYOL1O2) &  
   + C(ind_NOXYQO2) + C(ind_OXYQO2) + C(ind_OXNNCATCO2) + &  
   C(ind_OXNCATECO2) + C(ind_C534O2) + C(ind_CHOMOHCO3) + C(ind_CO2H3CO3) &  
   + C(ind_NMXYOL1O2) + C(ind_NMXYQO2) + C(ind_MXYQO2) + &  
   C(ind_MXNNCATCO2) + C(ind_MXNCATECO2) + C(ind_NPXYOL1O2) + &  
   C(ind_NPXYQO2) + C(ind_PXYQO2) + C(ind_PXNNCATCO2) + C(ind_PXNCATECO2) &  
   + C(ind_NPHEN1O2) + C(ind_NNCATECO2) + C(ind_NCATECO2) + C(ind_NBZQO2) &  
   + C(ind_PBZQO2) + C(ind_NPTLQO2) + C(ind_PTLQO2) + C(ind_NCRES1O2) + &  
   C(ind_MNNCATECO2) + C(ind_MNCATECO2) + C(ind_HOCH2CO3) + &  
   C(ind_CH3COCH2O2) + C(ind_MEKBO2) + C(ind_C6O4KETO2) + &  
   C(ind_NDNOXYOLO2) + C(ind_DNOXYOLO2) + C(ind_ACCOCOMEO2) + &  
   C(ind_C535O2) + C(ind_MMALANHYO2) + C(ind_CH3COCO3) + C(ind_MALANHYO2) &  
   + C(ind_NDNMXYOLO2) + C(ind_DNMXYOLO2) + C(ind_TL4OHNO2O2) + &  
   C(ind_NDNPXYOLO2) + C(ind_DNPXYOLO2) + C(ind_NDNPHENO2) + &  
   C(ind_DNPHENO2) + C(ind_NDNCRESO2) + C(ind_DNCRESO2) + &  
   C(ind_ACCOMECO3) + C(ind_C5M2OHOCO3) + C(ind_C4COMOHCO3) + &  
   C(ind_C23O3MCO3) + C(ind_C23O3CCO3) + C(ind_C7CO2OHCO3) + &  
   C(ind_C6MOHCOCO3) + C(ind_C5CO2OHCO3) + C(ind_C6CO2OHCO3) + &  
   C(ind_C4DBM2CO3) + C(ind_C5CO234O2) + C(ind_C5CO2DBCO3) + &  
   C(ind_C5DBCO2CO3) + C(ind_C4CO2DBCO3 ) + C(ind_MO2m) 


! End INLINED RCONST

  RCONST(1) = (7.45E4*(TEMP/300.)**-2.6)
  RCONST(2) = (8.0E-12*EXP(-2060./TEMP))
  RCONST(3) = (PAN(1.0E-31,0.0,-1.6,5.0E-11,0.0,-0.3,0.85))
  RCONST(4) = (5.5E-12*EXP(188./TEMP))
  RCONST(5) = (PAN(1.3E-31,0.0,-1.5,2.3E-11,0.0,0.24,0.6))
  RCONST(6) = (1.68E8*EXP(67./TEMP)+3.95E8*EXP(130./TEMP))
  RCONST(7) = (1.4E-12*EXP(-1310./TEMP))
  RCONST(8) = (1.4E-13*EXP(-2470./TEMP))
  RCONST(9) = (1.7E-20*EXP(530./TEMP))
  RCONST(10) = (1.8E-11*EXP(110./TEMP))
  RCONST(11) = (4.50E-14*EXP(-1260./TEMP))
  RCONST(12) = (PAN(3.6E-30,0.0,-4.1,1.9E-12,0.0,0.2,0.35))
! RCONST(13) = constant rate coefficient
  RCONST(14) = (1.70E-12*EXP(-940./TEMP))
  RCONST(15) = (7.7E-12*EXP(-2100./TEMP))
  RCONST(16) = (2.9E-12*EXP(-160./TEMP))
  RCONST(17) = (2.03E-16*EXP(693./TEMP)*(TEMP/300.)**4.57)
  RCONST(18) = (4.8E-11*EXP(250./TEMP))
  RCONST(19) = ((1+(3.5E-4*EXP(2200./TEMP)))*(2.20E-13*EXP(600./TEMP)+4.75E-14*EXP(980./TEMP)))
  RCONST(20) = (PAN(7.4E-31,0.0,-2.4,3.3E-11,0.0,-0.3,0.81))
  RCONST(21) = (PAN(3.2E-30,0.0,-4.5,3.0E-11,0.0,0.0,0.41))
! RCONST(22) = constant rate coefficient
  RCONST(23) = (3.45E-12*EXP(270./TEMP))
  RCONST(24) = (PAN(1.4E-31,0.0,-3.1,4.0E-12,0.0,0.0,0.4))
  RCONST(25) = (3.2E-13*EXP(690./TEMP))
! RCONST(26) = constant rate coefficient
  RCONST(27) = (2.5E-12*EXP(260./TEMP))
  RCONST(28) = (2.4E-14*EXP(460./TEMP)+(1.625E-14*EXP(1335./TEMP))/(1+(602.*EXP(-864./TEMP))))
  RCONST(29) = (SUN*3.78E-5)
  RCONST(30) = (SUN*4.41E-4)
  RCONST(31) = (SUN*7.88E-6)
  RCONST(32) = (SUN*8.92E-3)
  RCONST(33) = (SUN*2.23E-2)
  RCONST(34) = (SUN*1.54E-1)
  RCONST(35) = (SUN*1.98E-3)
  RCONST(36) = (SUN*6.85E-7)
  RCONST(37) = (GCARR(8.50E-13,0.0E+00,-2450.0))
  RCONST(38) = (SUN*4.5E-5)
  RCONST(39) = (SUN*5.21E-6)
  RCONST(40) = (PAN(1.3E-3,-1.1E4,-3.5,9.7E14,-1.108E4,0.1,0.35))
  RCONST(41) = (PAN(4.1E-5,-10650.0,0.0,6.0E15,-11170.0,0.0,0.4))
  RCONST(42) = (GCARR(2.80E-12,0.0E+00,300.0))
  RCONST(43) = (GCARR(4.10E-13,0.0E+00,750.0))
  RCONST(44) = (GC_TBRANCH(9.50E-14,0.0E+00,390.0,2.62E1,0.0,-1130.0))
  RCONST(45) = (GC_TBRANCH(9.50E-14,0.0E+00,390.0,4.E-2,0.0,1130.0))
  RCONST(46) = (GCJPLPR(1.00E-30,4.8E+00,0.0,7.2E-12,2.1E0,0.0,0.6,0.0,0.0))
  RCONST(47) = (GCJPLPR(1.05E-02,4.8E+00,-11234.0,7.58E16,2.1E0,-11234.0,0.6,0.0,0.0))
  RCONST(48) = (GCARR(2.66E-12,0.0E+00,200.0))
  RCONST(49) = (GCARR(1.14E-12,0.0E+00,200.0))
  RCONST(50) = (GCARR(5.50E-12,0.0E+00,125.0))
  RCONST(51) = (GCARR(5.80E-16,0.0E+00,0.0))
  RCONST(52) = (GCARR(2.90E-12,0.0E+00,-345.0))
  RCONST(53) = (KCO(1.44E-13,4.2E19))
  RCONST(54) = (SUN*6.49E-6)
  RCONST(55) = (SUN*4.65E-5)
  RCONST(56) = (SUN*6.01E-5)
  RCONST(57) = (SUN*0.05*1.95E-4)
  RCONST(58) = (SUN*0.95*1.95E-4)
  RCONST(59) = (GCARR(2.45E-12,0.0E+00,-1775.0))
  RCONST(60) = (GCARR(2.80E-12,0.0E+00,300.0))
  RCONST(61) = (GCARR(4.10E-13,0.0E+00,750.0))
  RCONST(62) = (GC_TBRANCH(9.50E-14,0.0E+00,390.0,2.62E1,0.0,-1130.0))
  RCONST(63) = (GC_TBRANCH(9.50E-14,0.0E+00,390.0,4.E-2,0.0,1130.0))
  RCONST(64) = (GC_TBRANCH(9.50E-14,0.0E+00,390.0,2.62E1,0.0,-1130.0))
  RCONST(65) = (GC_TBRANCH(9.50E-14,0.0E+00,390.0,4.E-2,0.0,1130.0))
  RCONST(66) = (GCJPLPR(1.00E-30,4.8E+00,0.0,7.2E-12,2.1E0,0.0,0.6,0.0,0.0))
  RCONST(67) = (GCJPLPR(1.05E-02,4.8E+00,-11234.0,7.58E16,2.1E0,-11234.0,0.6,0.0,0.0))
  RCONST(68) = (GCARR(2.66E-12,0.0E+00,200.0))
  RCONST(69) = (GCARR(1.14E-12,0.0E+00,200.0))
  RCONST(70) = (GCARR(5.50E-12,0.0E+00,125.0))
  RCONST(71) = (GCARR(5.80E-16,0.0E+00,0.0))
  RCONST(72) = (GCARR(2.90E-12,0.0E+00,-345.0))
  RCONST(73) = (KCO(1.44E-13,4.2E19))
  RCONST(74) = (SUN*6.49E-6)
  RCONST(75) = (SUN*4.65E-5)
  RCONST(76) = (SUN*6.01E-5)
  RCONST(77) = (SUN*0.05*1.95E-4)
  RCONST(78) = (SUN*0.95*1.95E-4)
  RCONST(79) = (SUN*3.1E-3)
  RCONST(80) = (GCARR(3.36E-12,0.0E+00,-1860.0))
  RCONST(81) = (GCARR(1.50E-11,0.0E+00,0.0))
  RCONST(82) = (SUN*8.46E-5)
  RCONST(83) = (SUN*8.61E-6)
  RCONST(84) = (SUN*3.71E-5)
  RCONST(85) = (GC_GLYXNO3(1.40E-12,0.0E+00,-1860.0))
  RCONST(86) = (GCARR(3.10E-12,0.0E+00,340.0))
  RCONST(87) = (SUN*6.4E-6)
  RCONST(88) = (GCARR(4.63E-12,0.0E+00,350.0))
  RCONST(89) = (GCARR(1.40E-12,0.0E+00,-1900.0))
  RCONST(90) = (GCARR(6.13E-13,0.0E+00,200.0))
  RCONST(91) = (SUN*3.0E-5)
  RCONST(92) = (GCJPLEQ(9.30E-29,0.0E+00,14000.0,9.7E-29,5.6E0,0.0,9.3E-12,1.5E0,0.,0.6,0.,0.))
  RCONST(93) = (SUN*1.0E-6)
  RCONST(94) = (GCARR(1.80E-12,0.0E+00,500.0)*0.1*RO2)
  RCONST(95) = (GCARR(1.80E-12,0.0E+00,500.0)*0.9*RO2)
  RCONST(96) = (GCJPLPR(9.70E-29,5.6E+00,0.0,9.3E-12,1.5E0,0.0,0.6,0.0,0.0))
  RCONST(97) = (GCARR(8.10E-12,0.0E+00,270.0))
  RCONST(98) = (GCARR(5.20E-13,0.0E+00,980.0))
  RCONST(99) = (GC_GLYCOHA(8.00E-12,0.0E+00,0.0))
  RCONST(100) = (GC_GLYCOHB(8.00E-12,0.0E+00,0.0))
  RCONST(101) = (SUN*1.67E-5)
  RCONST(102) = (4.10D-16*0.33)
  RCONST(103) = (1.36D-11*0.55*0.33)
  RCONST(104) = (1.36D-11*0.24*0.33)
  RCONST(105) = (1.36D-11*0.05*0.33)
  RCONST(106) = (1.36D-11*0.16*0.33)
  RCONST(107) = (2.60D-16*0.33)
  RCONST(108) = (2.31D-11*0.50*0.33)
  RCONST(109) = (2.31D-11*0.29*0.33)
  RCONST(110) = (2.31D-11*0.04*0.33)
  RCONST(111) = (2.31D-11*0.17*0.33)
  RCONST(112) = (5.00D-16*0.34)
  RCONST(113) = (1.43D-11*0.625*0.34)
  RCONST(114) = (1.43D-11*0.155*0.34)
  RCONST(115) = (1.43D-11*0.10*0.34)
  RCONST(116) = (1.43D-11*0.12*0.34)
  RCONST(117) = (2.3D-12*EXP(-190/TEMP)*0.352)
  RCONST(118) = (2.3D-12*EXP(-190/TEMP)*0.118)
  RCONST(119) = (2.3D-12*EXP(-190/TEMP)*0.53)
  RCONST(120) = (1.8D-12*EXP(340/TEMP)*0.07)
  RCONST(121) = (1.8D-12*EXP(340/TEMP)*0.18)
  RCONST(122) = (1.8D-12*EXP(340/TEMP)*0.65)
  RCONST(123) = (1.8D-12*EXP(340/TEMP)*0.10)
  RCONST(124) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(125) = (2.7D-12*EXP(360/TEMP)*0.105)
  RCONST(126) = (2.7D-12*EXP(360/TEMP)*0.895)
! RCONST(127) = constant rate coefficient
  RCONST(128) = (1.30D-12*0.20*RO2)
  RCONST(129) = (1.30D-12*0.60*RO2)
  RCONST(130) = (1.30D-12*0.20*RO2)
  RCONST(131) = (2.7D-12*EXP(360/TEMP)*0.138)
  RCONST(132) = (2.7D-12*EXP(360/TEMP)*0.862)
! RCONST(133) = constant rate coefficient
  RCONST(134) = (9.20D-14*RO2*0.70)
  RCONST(135) = (9.20D-14*RO2*0.30)
  RCONST(136) = (2*1.4D-12*EXP(-1860/TEMP)*4.25)
  RCONST(137) = (2.00D-18)
  RCONST(138) = (5.99D-11*0.632)
  RCONST(139) = (5.99D-11*0.368)
  RCONST(140) = (SUN*8.92E-3*0.05*0.5)
  RCONST(141) = (SUN*8.92E-3*0.05*0.5)
  RCONST(142) = (3.20D-11*0.51)
  RCONST(143) = (3.20D-11*0.39)
  RCONST(144) = (3.20D-11*0.10)
  RCONST(145) = (8.00D-11*0.07)
  RCONST(146) = (8.00D-11*0.73)
  RCONST(147) = (8.00D-11*0.20)
  RCONST(148) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(149) = (2.7D-12*EXP(360/TEMP)*0.105)
  RCONST(150) = (2.7D-12*EXP(360/TEMP)*0.895)
! RCONST(151) = constant rate coefficient
  RCONST(152) = (1.30D-12*0.20*RO2)
  RCONST(153) = (1.30D-12*0.60*RO2)
  RCONST(154) = (1.30D-12*0.20*RO2)
  RCONST(155) = (2.7D-12*EXP(360/TEMP)*0.138)
  RCONST(156) = (2.7D-12*EXP(360/TEMP)*0.862)
! RCONST(157) = constant rate coefficient
  RCONST(158) = (8.80D-13*RO2*0.6)
  RCONST(159) = (8.80D-13*RO2*0.2)
  RCONST(160) = (8.80D-13*RO2*0.2)
  RCONST(161) = (1.4D-12*EXP(-1860/TEMP)*4.25)
  RCONST(162) = (2.00D-18)
  RCONST(163) = (8.02D-11*0.26)
  RCONST(164) = (8.02D-11*0.74)
  RCONST(165) = (SUN*8.92E-3*0.028*0.5)
  RCONST(166) = (SUN*8.92E-3*0.028*0.5)
  RCONST(167) = (4.50D-11*0.39)
  RCONST(168) = (4.50D-11*0.10)
  RCONST(169) = (4.50D-11*0.51)
  RCONST(170) = (9.10D-11*0.07)
  RCONST(171) = (9.10D-11*0.42)
  RCONST(172) = (9.10D-11*0.51)
  RCONST(173) = (2.91D-13*EXP(1300/TEMP)*0.820)
  RCONST(174) = (2.7D-12*EXP(360/TEMP)*0.105)
  RCONST(175) = (2.7D-12*EXP(360/TEMP)*0.895)
! RCONST(176) = constant rate coefficient
  RCONST(177) = (1.30D-12*0.20*RO2)
  RCONST(178) = (1.30D-12*0.60*RO2)
  RCONST(179) = (1.30D-12*0.20*RO2)
  RCONST(180) = (2.7D-12*EXP(360/TEMP)*0.138)
  RCONST(181) = (2.7D-12*EXP(360/TEMP)*0.862)
! RCONST(182) = constant rate coefficient
  RCONST(183) = (9.20D-14*RO2*0.70)
  RCONST(184) = (9.20D-14*RO2*0.30)
  RCONST(185) = (1.4D-12*EXP(-1860/TEMP)*4.25)
  RCONST(186) = (5.00D-18)
  RCONST(187) = (7.88D-11*0.26)
  RCONST(188) = (7.88D-11*0.74)
  RCONST(189) = (SUN*8.92E-3*0.1*0.5)
  RCONST(190) = (SUN*8.92E-3*0.1*0.5)
  RCONST(191) = (3.48D-11*0.51)
  RCONST(192) = (3.48D-11*0.39)
  RCONST(193) = (3.48D-11*0.10)
  RCONST(194) = (8.00D-11*0.07)
  RCONST(195) = (8.00D-11*0.73)
  RCONST(196) = (8.00D-11*0.2)
  RCONST(197) = (2.7D-12*EXP(360/TEMP)*0.082)
  RCONST(198) = (2.7D-12*EXP(360/TEMP)*0.918)
! RCONST(199) = constant rate coefficient
  RCONST(200) = (8.80D-13*RO2*0.2)
  RCONST(201) = (8.80D-13*RO2*0.6)
  RCONST(202) = (8.80D-13*RO2*0.2)
  RCONST(203) = (2*1.4D-12*EXP(-1860/TEMP)*2.75)
  RCONST(204) = (2.00D-18)
  RCONST(205) = (6.08D-11*0.31)
  RCONST(206) = (6.08D-11*0.69)
  RCONST(207) = (SUN*8.92E-3*0.1*0.5)
  RCONST(208) = (SUN*8.92E-3*0.1*0.5)
  RCONST(209) = (3.8D-12*0.742)
  RCONST(210) = (3.8D-12*0.258)
  RCONST(211) = (4.7D-13*EXP(1220/TEMP)*0.06)
  RCONST(212) = (4.7D-13*EXP(1220/TEMP)*0.8)
  RCONST(213) = (4.7D-13*EXP(1220/TEMP)*0.14)
  RCONST(214) = (1.5D-13*EXP(1310/TEMP))
  RCONST(215) = (2.7D-12*EXP(360/TEMP)*0.105)
  RCONST(216) = (2.7D-12*EXP(360/TEMP)*0.895)
! RCONST(217) = constant rate coefficient
  RCONST(218) = (2*(1.03D-13*EXP(365/TEMP)*2.4D-14*EXP(1620/TEMP))**0.5*RO2*0.2)
  RCONST(219) = (2*(1.03D-13*EXP(365/TEMP)*2.4D-14*EXP(1620/TEMP))**0.5*RO2*0.6)
  RCONST(220) = (2*(1.03D-13*EXP(365/TEMP)*2.4D-14*EXP(1620/TEMP))**0.5*RO2*0.2)
  RCONST(221) = (1.4D-11*0.103)
  RCONST(222) = (1.4D-11*0.506)
  RCONST(223) = (1.4D-11*0.391)
  RCONST(224) = (4.65D-11*0.200)
  RCONST(225) = (4.65D-11*0.727)
  RCONST(226) = (4.65D-11*0.073)
  RCONST(227) = (2.7D-12*EXP(360/TEMP)*0.889)
  RCONST(228) = (2.7D-12*EXP(360/TEMP)*0.111)
! RCONST(229) = constant rate coefficient
  RCONST(230) = (8.80D-13*RO2*0.20)
  RCONST(231) = (8.80D-13*RO2*0.60)
  RCONST(232) = (8.80D-13*RO2*0.20)
  RCONST(233) = (1.4D-12*EXP(-1860/TEMP)*2.75)
  RCONST(234) = (5.00D-18)
  RCONST(235) = (7.99D-11*0.31)
  RCONST(236) = (7.99D-11*0.69)
  RCONST(237) = (SUN*8.92E-3*0.1*0.5)
  RCONST(238) = (SUN*8.92E-3*0.1*0.5)
  RCONST(239) = (2.75D-11)
  RCONST(240) = (SUN*5.79E-6)
  RCONST(241) = (1.31D-11)
  RCONST(242) = (SUN*1.79E-6)
  RCONST(243) = (2.50D-14*EXP(-300/TEMP))
  RCONST(244) = (1.4D-12*EXP(-1860/TEMP))
  RCONST(245) = (1.29D-11)
  RCONST(246) = (SUN*1.10E-6)
  RCONST(247) = (SUN*1.10E-6)
  RCONST(248) = (1.60D-11)
  RCONST(249) = (7.62D-11)
  RCONST(250) = (SUN*5.79E-6)
  RCONST(251) = (7.28D-11)
  RCONST(252) = (SUN*8.33E-6)
  RCONST(253) = (1.00D+06*0.175)
  RCONST(254) = (1.00D+06*0.245)
  RCONST(255) = (1.00D+06*0.175)
  RCONST(256) = (1.00D+06*0.08)
  RCONST(257) = (1.00D+06*0.08)
  RCONST(258) = (1.00D+06*0.245)
  RCONST(259) = (1.11D-10)
  RCONST(260) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(261) = (5.2D-13*EXP(980/TEMP)*0.15)
  RCONST(262) = (5.2D-13*EXP(980/TEMP)*0.41)
  RCONST(263) = (7.5D-12*EXP(290/TEMP))
  RCONST(264) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(265) = (2.3E-12*1.74)
  RCONST(266) = (1.00D-11*RO2*0.70)
  RCONST(267) = (1.00D-11*RO2*0.30)
  RCONST(268) = (2*1.4D-12*EXP(-1860/TEMP)*8.5)
  RCONST(269) = (4.24D-11)
  RCONST(270) = (SUN*5.50E-5*2)
  RCONST(271) = (1.00D+06*0.18)
  RCONST(272) = (1.00D+06*0.125)
  RCONST(273) = (1.00D+06*0.125)
  RCONST(274) = (1.00D+06*0.57)
  RCONST(275) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(276) = (2.7D-12*EXP(360/TEMP)*0.173)
  RCONST(277) = (2.7D-12*EXP(360/TEMP)*0.827)
! RCONST(278) = constant rate coefficient
  RCONST(279) = (8.80D-12*RO2*0.20)
  RCONST(280) = (8.80D-12*RO2*0.60)
  RCONST(281) = (8.80D-12*RO2*0.20)
  RCONST(282) = (2.91D-13*EXP(1300/TEMP)*0.820)
  RCONST(283) = (2.7D-12*EXP(360/TEMP))
! RCONST(284) = constant rate coefficient
  RCONST(285) = (9.20D-14*RO2*0.70)
  RCONST(286) = (9.20D-14*RO2*0.30)
  RCONST(287) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(288) = (2.7D-12*EXP(360/TEMP))
! RCONST(289) = constant rate coefficient
  RCONST(290) = (8.00D-13*RO2*0.70)
  RCONST(291) = (8.00D-13*RO2*0.30)
  RCONST(292) = (2.08D-12)
  RCONST(293) = (2.86D-13)
  RCONST(294) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(295) = (2.7D-12*EXP(360/TEMP))
! RCONST(296) = constant rate coefficient
  RCONST(297) = (8.00D-13*RO2*0.70)
  RCONST(298) = (8.00D-13*RO2*0.30)
  RCONST(299) = (2.01D-10)
  RCONST(300) = (2.81D-17)
  RCONST(301) = (2.05D-10)
  RCONST(302) = (2.05D-11)
  RCONST(303) = (SUN*5.79E-6)
  RCONST(304) = (6.03D-12)
  RCONST(305) = (SUN*1.79E-6)
  RCONST(306) = (2.50D-14*EXP(-300/TEMP))
  RCONST(307) = (1.4D-12*EXP(-1860/TEMP))
  RCONST(308) = (1.29D-11)
  RCONST(309) = (SUN*1.10E-6)
  RCONST(310) = (SUN*1.10E-6)
  RCONST(311) = (8.93D-12)
  RCONST(312) = (1.17D-10)
  RCONST(313) = (SUN*5.79E-6)
  RCONST(314) = (7.15D-11)
  RCONST(315) = (SUN*2.99E-6)
  RCONST(316) = (1.00D+06*0.11)
  RCONST(317) = (1.00D+06*0.11)
  RCONST(318) = (1.00D+06*0.26)
  RCONST(319) = (1.00D+06*0.26)
  RCONST(320) = (1.00D+06*0.26)
  RCONST(321) = (2.50D-10)
  RCONST(322) = (7.93D-11*0.5)
  RCONST(323) = (7.93D-11*0.5)
  RCONST(324) = (SUN*3.98E-6*0.5)
  RCONST(325) = (SUN*3.98E-6*0.5)
  RCONST(326) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(327) = (5.2D-13*EXP(980/TEMP)*0.15)
  RCONST(328) = (5.2D-13*EXP(980/TEMP)*0.41)
  RCONST(329) = (7.5D-12*EXP(290/TEMP))
  RCONST(330) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(331) = (2.3E-12*1.74)
  RCONST(332) = (1.00D-11*RO2*0.70)
  RCONST(333) = (1.00D-11*RO2*0.30)
  RCONST(334) = (2*1.4D-12*EXP(-1860/TEMP)*8.5)
  RCONST(335) = (4.28D-11)
  RCONST(336) = (SUN*5.50E-5*2)
  RCONST(337) = (1.00D+06*0.125)
  RCONST(338) = (1.00D+06*0.57)
  RCONST(339) = (1.00D+06*0.125)
  RCONST(340) = (1.00D+06*0.18)
  RCONST(341) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(342) = (2.7D-12*EXP(360/TEMP)*0.173)
  RCONST(343) = (2.7D-12*EXP(360/TEMP)*0.827)
! RCONST(344) = constant rate coefficient
  RCONST(345) = (8.80D-13*RO2*0.2)
  RCONST(346) = (8.80D-13*RO2*0.6)
  RCONST(347) = (8.80D-13*RO2*0.2)
  RCONST(348) = (2.91D-13*EXP(1300/TEMP)*0.820)
  RCONST(349) = (2.7D-12*EXP(360/TEMP))
! RCONST(350) = constant rate coefficient
  RCONST(351) = (8.80D-13*RO2*0.2)
  RCONST(352) = (8.80D-13*RO2*0.2)
  RCONST(353) = (8.80D-13*RO2*0.6)
  RCONST(354) = (2.08D-12)
  RCONST(355) = (2.86D-13)
  RCONST(356) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(357) = (2.7D-12*EXP(360/TEMP))
! RCONST(358) = constant rate coefficient
  RCONST(359) = (8.00D-13*RO2*0.70)
  RCONST(360) = (8.00D-13*RO2*0.30)
  RCONST(361) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(362) = (2.7D-12*EXP(360/TEMP))
! RCONST(363) = constant rate coefficient
  RCONST(364) = (8.00D-13*RO2*0.70)
  RCONST(365) = (8.00D-13*RO2*0.30)
  RCONST(366) = (2.01D-10)
  RCONST(367) = (2.81D-17)
  RCONST(368) = (2.05D-10)
  RCONST(369) = (2.05D-11)
  RCONST(370) = (SUN*5.79E-6)
  RCONST(371) = (6.03D-12)
  RCONST(372) = (SUN*1.79E-6)
  RCONST(373) = (2.50D-14*EXP(-300/TEMP))
  RCONST(374) = (1.4D-12*EXP(-1860/TEMP))
  RCONST(375) = (1.29D-11)
  RCONST(376) = (SUN*1.10E-6)
  RCONST(377) = (SUN*1.10E-6)
  RCONST(378) = (8.93D-12)
  RCONST(379) = (7.65D-11)
  RCONST(380) = (SUN*5.79E-6)
  RCONST(381) = (7.16D-11)
  RCONST(382) = (SUN*8.33E-6)
  RCONST(383) = (1.00D+06*0.64)
  RCONST(384) = (1.00D+06*0.18)
  RCONST(385) = (1.00D+06*0.18)
  RCONST(386) = (1.12D-10)
  RCONST(387) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(388) = (5.2D-13*EXP(980/TEMP)*0.15)
  RCONST(389) = (5.2D-13*EXP(980/TEMP)*0.41)
  RCONST(390) = (7.5D-12*EXP(290/TEMP))
  RCONST(391) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(392) = (2.3E-12*1.74)
  RCONST(393) = (1.00D-11*RO2*0.70)
  RCONST(394) = (1.00D-11*RO2*0.30)
  RCONST(395) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(396) = (2.7D-12*EXP(360/TEMP)*0.173)
  RCONST(397) = (2.7D-12*EXP(360/TEMP)*0.827)
! RCONST(398) = constant rate coefficient
  RCONST(399) = (8.80D-12*RO2*0.2)
  RCONST(400) = (8.80D-12*RO2*0.6)
  RCONST(401) = (8.80D-12*RO2*0.2)
  RCONST(402) = (2.91D-13*EXP(1300/TEMP)*0.820)
  RCONST(403) = (2.7D-12*EXP(360/TEMP))
! RCONST(404) = constant rate coefficient
  RCONST(405) = (9.20D-14*RO2*0.70)
  RCONST(406) = (9.20D-14*RO2*0.30)
  RCONST(407) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(408) = (2.7D-12*EXP(360/TEMP))
! RCONST(409) = constant rate coefficient
  RCONST(410) = (8.00D-13*RO2*0.70)
  RCONST(411) = (8.00D-13*RO2*0.30)
  RCONST(412) = (2.08D-12)
  RCONST(413) = (2.86D-13)
  RCONST(414) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(415) = (2.7D-12*EXP(360/TEMP))
! RCONST(416) = constant rate coefficient
  RCONST(417) = (8.00D-13*RO2*0.70)
  RCONST(418) = (8.00D-13*RO2*0.30)
  RCONST(419) = (2.01D-10)
  RCONST(420) = (2.81D-17)
  RCONST(421) = (2.05D-10)
  RCONST(422) = (9.77D-11)
  RCONST(423) = (SUN*5.79E-6)
  RCONST(424) = (7.30D-11)
  RCONST(425) = (SUN*2.99E-6)
  RCONST(426) = (1.00D+06*0.50)
  RCONST(427) = (1.00D+06*0.50)
  RCONST(428) = (1.21D-10)
  RCONST(429) = (8.16D-11)
  RCONST(430) = (SUN*3.98E-6)
  RCONST(431) = (5.2D-13*EXP(980/TEMP)*0.15)
  RCONST(432) = (5.2D-13*EXP(980/TEMP)*0.41)
  RCONST(433) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(434) = (7.5D-12*EXP(290/TEMP))
  RCONST(435) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(436) = (2.3E-12*1.74)
  RCONST(437) = (1.00D-11*RO2*0.3)
  RCONST(438) = (1.00D-11*RO2*0.7)
  RCONST(439) = (2*1.4D-12*EXP(-1860/TEMP)*4.0)
  RCONST(440) = (4.32D-11)
  RCONST(441) = (SUN*5.50E-5*2)
  RCONST(442) = (2.91D-13*EXP(1300/TEMP)*0.77)
  RCONST(443) = (2.7D-12*EXP(360/TEMP)*0.105)
  RCONST(444) = (2.7D-12*EXP(360/TEMP)*0.895)
! RCONST(445) = constant rate coefficient
  RCONST(446) = (8.80D-13*0.2*RO2)
  RCONST(447) = (8.80D-13*0.6*RO2)
  RCONST(448) = (8.80D-13*0.2*RO2)
  RCONST(449) = (2.91D-13*EXP(1300/TEMP)*0.706)
  RCONST(450) = (2.7D-12*EXP(360/TEMP))
! RCONST(451) = constant rate coefficient
  RCONST(452) = (8.80D-13*RO2*0.2)
  RCONST(453) = (8.80D-13*RO2*0.6)
  RCONST(454) = (8.80D-13*RO2*0.2)
  RCONST(455) = (2*1.4D-12*EXP(-1860/TEMP)*2.0)
  RCONST(456) = (2.00D-18)
  RCONST(457) = (5.20D-11*0.83)
  RCONST(458) = (5.20D-11*0.17)
  RCONST(459) = (SUN*8.92E-3*0.14*0.4)
  RCONST(460) = (SUN*8.92E-3*0.14*0.6)
  RCONST(461) = (2.08D-12)
  RCONST(462) = (2.86D-13)
  RCONST(463) = (2.91D-13*EXP(1300/TEMP)*0.770)
  RCONST(464) = (2.7D-12*EXP(360/TEMP))
! RCONST(465) = constant rate coefficient
  RCONST(466) = (8.00D-13*RO2*0.7)
  RCONST(467) = (8.00D-13*RO2*0.3)
  RCONST(468) = (9.9D-11)
  RCONST(469) = (9.2D-18)
  RCONST(470) = (1.0D-10)
  RCONST(471) = (2.91D-13*EXP(1300/TEMP)*0.770)
  RCONST(472) = (2.7D-12*EXP(360/TEMP))
! RCONST(473) = constant rate coefficient
  RCONST(474) = (8.00D-13*RO2*0.7)
  RCONST(475) = (8.00D-13*RO2*0.3)
  RCONST(476) = (2.05D-11)
  RCONST(477) = (SUN*5.79E-6)
  RCONST(478) = (6.03D-12)
  RCONST(479) = (SUN*1.79E-6)
  RCONST(480) = (2.50D-14*EXP(-300/TEMP))
  RCONST(481) = (SUN*3.30E-6)
  RCONST(482) = (2.40D-15)
  RCONST(483) = (5.9D-12*EXP(225/TEMP))
  RCONST(484) = (2.7D-11)
  RCONST(485) = (2.91D-13*EXP(1300/TEMP)*0.820)
  RCONST(486) = (2.7D-12*EXP(360/TEMP))
! RCONST(487) = constant rate coefficient
  RCONST(488) = (8.00D-13*RO2*0.70)
  RCONST(489) = (8.00D-13*RO2*0.30)
  RCONST(490) = (2.91D-13*EXP(1300/TEMP)*0.820)
  RCONST(491) = (2.7D-12*EXP(360/TEMP))
! RCONST(492) = constant rate coefficient
  RCONST(493) = (8.00D-13*RO2*0.70)
  RCONST(494) = (8.00D-13*RO2*0.30)
  RCONST(495) = (2.08D-12)
  RCONST(496) = (2.86D-13)
  RCONST(497) = (1.7D-10*1.0)
  RCONST(498) = (2.8D-17)
  RCONST(499) = (2.0D-10*1.0)
  RCONST(500) = (9.64D-11)
  RCONST(501) = (SUN*5.79E-6)
  RCONST(502) = (1.00D+06*0.2)
  RCONST(503) = (1.00D+06*0.2)
  RCONST(504) = (1.00D+06*0.2)
  RCONST(505) = (1.00D+06*0.2)
  RCONST(506) = (1.00D+06*0.2)
  RCONST(507) = (7.16D-11)
  RCONST(508) = (SUN*2.99E-6)
  RCONST(509) = (1.20D-10)
  RCONST(510) = (7.99D-11)
  RCONST(511) = (SUN*3.98E-6)
  RCONST(512) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(513) = (5.2D-13*EXP(980/TEMP)*0.15)
  RCONST(514) = (5.2D-13*EXP(980/TEMP)*0.41)
  RCONST(515) = (7.5D-12*EXP(290/TEMP))
  RCONST(516) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(517) = (2.3E-12*1.74)
  RCONST(518) = (1.00D-11*RO2*0.7)
  RCONST(519) = (1.00D-11*RO2*0.3)
  RCONST(520) = (2.91D-13*EXP(1300/TEMP)*0.820)
  RCONST(521) = (2.7D-12*EXP(360/TEMP)*0.105)
  RCONST(522) = (2.7D-12*EXP(360/TEMP)*0.895)
! RCONST(523) = constant rate coefficient
  RCONST(524) = (8.80D-13*RO2*0.20)
  RCONST(525) = (8.80D-13*RO2*0.60)
  RCONST(526) = (8.80D-13*RO2*0.20)
  RCONST(527) = (2.91D-13*EXP(1300/TEMP)*0.770)
  RCONST(528) = (2.7D-12*EXP(360/TEMP))
! RCONST(529) = constant rate coefficient
  RCONST(530) = (8.80D-13*RO2*0.2)
  RCONST(531) = (8.80D-13*RO2*0.6)
  RCONST(532) = (8.80D-13*RO2*0.2)
  RCONST(533) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(534) = (5.2D-13*EXP(980/TEMP)*0.15)
  RCONST(535) = (5.2D-13*EXP(980/TEMP)*0.41)
  RCONST(536) = (7.5D-12*EXP(290/TEMP))
  RCONST(537) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(538) = (2.3E-12*1.74)
  RCONST(539) = (1.00D-11*0.70*RO2)
  RCONST(540) = (1.00D-11*0.30*RO2)
  RCONST(541) = (2.91D-13*EXP(1300/TEMP)*0.820)
  RCONST(542) = (2.7D-12*EXP(360/TEMP))
! RCONST(543) = constant rate coefficient
  RCONST(544) = (2.50D-13*RO2)
  RCONST(545) = (8.06D-11)
  RCONST(546) = (SUN*3.98E-6)
  RCONST(547) = (3.00D-13)
  RCONST(548) = (2.20D-19)
  RCONST(549) = (4.45D-11)
  RCONST(550) = (SUN*2.68E-4)
  RCONST(551) = (1.40D-18*TEMP**2*EXP(194/TEMP))
  RCONST(552) = (1.4D-12*EXP(-1860/TEMP)*2.75)
  RCONST(553) = (2.00D-18*0.5)
  RCONST(554) = (2.00D-18*0.5)
  RCONST(555) = (6.2D-11*0.48)
  RCONST(556) = (6.2D-11*0.52)
  RCONST(557) = (SUN*8.92E-3*0.2*0.6)
  RCONST(558) = (SUN*8.92E-3*0.2*0.4)
  RCONST(559) = (1.4D-12*EXP(-1860/TEMP)*4.25)
  RCONST(560) = (5.00D-18*0.5)
  RCONST(561) = (5.00D-18*0.5)
  RCONST(562) = (7.10D-11*0.28)
  RCONST(563) = (7.10D-11*0.72)
  RCONST(564) = (SUN*8.92E-3*0.2*0.6)
  RCONST(565) = (SUN*8.92E-3*0.2*0.4)
  RCONST(566) = (1.00D-12)
  RCONST(567) = (1.40D-18)
  RCONST(568) = (3.01D-11)
  RCONST(569) = (1.00D-12)
  RCONST(570) = (8.00D-19)
  RCONST(571) = (6.90D-11)
  RCONST(572) = (3.98D-11)
  RCONST(573) = (SUN*1.10E-6+SUN*1.10E-6)
  RCONST(574) = (4.28D-11)
  RCONST(575) = (SUN*5.79E-6+SUN*1.10E-6+SUN*1.10E-6)
  RCONST(576) = (3.96D-11)
  RCONST(577) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(578) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(579) = (5.2D-13*EXP(980/TEMP)*0.15)
  RCONST(580) = (5.2D-13*EXP(980/TEMP)*0.41)
  RCONST(581) = (7.5D-12*EXP(290/TEMP))
  RCONST(582) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(583) = (2.3E-12*1.74)
  RCONST(584) = (1.00D-11*RO2*0.70)
  RCONST(585) = (1.00D-11*RO2*0.30)
  RCONST(586) = (2.91D-13*EXP(1300/TEMP)*0.706)
  RCONST(587) = (2.7D-12*EXP(360/TEMP))
! RCONST(588) = constant rate coefficient
  RCONST(589) = (9.20D-14*RO2)
  RCONST(590) = (1.2D-15)
  RCONST(591) = (1.0D-14)
  RCONST(592) = (1.0D-15)
  RCONST(593) = (7.0D-14)
  RCONST(594) = (6.0D-18)
  RCONST(595) = (1.0D-17)
  RCONST(596) = (9.33D-11)
  RCONST(597) = (SUN*5.79E-6+SUN*5.50E-5*2)
  RCONST(598) = (3.27D-11)
  RCONST(599) = (SUN*5.50E-5*2)
  RCONST(600) = (1.00D+06*0.5)
  RCONST(601) = (1.00D+06*0.5)
  RCONST(602) = (4.84D-11)
  RCONST(603) = (SUN*5.50E-5+SUN*1.25E-4)
  RCONST(604) = (8.11D-11)
  RCONST(605) = (SUN*5.50E-5*2)
  RCONST(606) = (4.28D-11)
  RCONST(607) = (SUN*1.10E-6+SUN*1.10E-6)
  RCONST(608) = (SUN*3.98E-6)
  RCONST(609) = (SUN*5.79E-6)
  RCONST(610) = (1.00D+06)
  RCONST(611) = (5.24D-11)
  RCONST(612) = (SUN*1.10E-6+SUN*1.10E-6)
  RCONST(613) = (SUN*3.98E-6)
  RCONST(614) = (1.43D-10)
  RCONST(615) = (SUN*5.79E-6)
  RCONST(616) = (SUN*2.99E-6)
  RCONST(617) = (1.00D+06)
  RCONST(618) = (1.40D-10)
  RCONST(619) = (SUN*2.99E-6)
  RCONST(620) = (3.13D-13)
  RCONST(621) = (1.53D-12)
  RCONST(622) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(623) = (2.7D-12*EXP(360/TEMP))
! RCONST(624) = constant rate coefficient
  RCONST(625) = (2.50D-13*RO2)
  RCONST(626) = (1.51D-10)
  RCONST(627) = (SUN*5.79E-6)
  RCONST(628) = (1.00D+06*0.52)
  RCONST(629) = (1.00D+06*0.48)
  RCONST(630) = (1.48D-10)
  RCONST(631) = (2.08D-12)
  RCONST(632) = (2.86D-13)
  RCONST(633) = (1.00D+06)
  RCONST(634) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(635) = (5.2D-13*EXP(980/TEMP)*0.15)
  RCONST(636) = (5.2D-13*EXP(980/TEMP)*0.41)
  RCONST(637) = (7.5D-12*EXP(290/TEMP))
  RCONST(638) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(639) = (2.3E-12*1.74)
  RCONST(640) = (1.00D-11*0.70*RO2)
  RCONST(641) = (1.00D-11*0.30*RO2)
  RCONST(642) = (2.91D-13*EXP(1300/TEMP)*0.820)
  RCONST(643) = (2.7D-12*EXP(360/TEMP))
! RCONST(644) = constant rate coefficient
  RCONST(645) = (2.50D-13*RO2)
  RCONST(646) = (1.4D-12*EXP(-1860/TEMP)*4.25)
  RCONST(647) = (2.0D-18*0.5)
  RCONST(648) = (2.0D-18*0.5)
  RCONST(649) = (7.00D-11*0.5)
  RCONST(650) = (7.00D-11*0.5)
  RCONST(651) = (SUN*8.92E-3*0.2*0.6)
  RCONST(652) = (SUN*8.92E-3*0.2*0.4)
  RCONST(653) = (1.00D-12)
  RCONST(654) = (8.00D-19)
  RCONST(655) = (2.67D-11)
  RCONST(656) = (1.4D-12*EXP(-1860/TEMP)*4.25)
  RCONST(657) = (1.4D-12*EXP(-1860/TEMP)*4.25)
  RCONST(658) = (5.00D-18*0.50)
  RCONST(659) = (5.00D-18*0.50)
  RCONST(660) = (4.41D-11*0.385)
  RCONST(661) = (4.41D-11*0.23)
  RCONST(662) = (4.41D-11*0.385)
  RCONST(663) = (SUN*8.92E-3*0.2*0.3)
  RCONST(664) = (SUN*8.92E-3*0.2*0.3)
  RCONST(665) = (SUN*8.92E-3*0.2*0.4)
  RCONST(666) = (1.00D-12)
  RCONST(667) = (8.00D-19)
  RCONST(668) = (2.42D-11)
  RCONST(669) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(670) = (5.2D-13*EXP(980/TEMP)*0.56)
  RCONST(671) = (7.5D-12*EXP(290/TEMP))
  RCONST(672) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(673) = (2.3E-12*1.74)
  RCONST(674) = (1.00D-11*RO2)
  RCONST(675) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(676) = (5.2D-13*EXP(980/TEMP)*0.15)
  RCONST(677) = (5.2D-13*EXP(980/TEMP)*0.41)
  RCONST(678) = (7.5D-12*EXP(290/TEMP))
  RCONST(679) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(680) = (2.3E-12*1.74)
  RCONST(681) = (1.00D-11*0.70*RO2)
  RCONST(682) = (1.00D-11*0.30*RO2)
  RCONST(683) = (6.00D-11)
  RCONST(684) = (SUN*1.80E-6+SUN*1.80E-6)
  RCONST(685) = (6.31D-11)
  RCONST(686) = (SUN*5.79E-6+SUN*1.80E-6+SUN*1.80E-6)
  RCONST(687) = (5.99D-11)
  RCONST(688) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(689) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(690) = (5.2D-13*EXP(980/TEMP)*0.15)
  RCONST(691) = (5.2D-13*EXP(980/TEMP)*0.41)
  RCONST(692) = (7.5D-12*EXP(290/TEMP))
  RCONST(693) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(694) = (2.3E-12*1.74)
  RCONST(695) = (1.00D-11*RO2*0.70)
  RCONST(696) = (1.00D-11*RO2*0.30)
  RCONST(697) = (2.91D-13*EXP(1300/TEMP)*0.625)
  RCONST(698) = (2.7D-12*EXP(360/TEMP))
! RCONST(699) = constant rate coefficient
  RCONST(700) = (9.20D-14*RO2*0.7)
  RCONST(701) = (9.20D-14*RO2*0.3)
  RCONST(702) = (1.2D-15)
  RCONST(703) = (1.0D-14)
  RCONST(704) = (1.0D-15)
  RCONST(705) = (7.0D-14)
  RCONST(706) = (1.0D-17)
  RCONST(707) = (6.0D-18)
  RCONST(708) = (7.41D-11)
  RCONST(709) = (SUN*5.79E-6+SUN*5.50E-5+SUN*3.98E-6)
  RCONST(710) = (3.08D-11)
  RCONST(711) = (SUN*5.50E-5+SUN*3.98E-6)
  RCONST(712) = (1.00D+06*0.5)
  RCONST(713) = (1.00D+06*0.5)
  RCONST(714) = (4.08D-11)
  RCONST(715) = (SUN*5.50E-5+SUN*2.68E-4)
  RCONST(716) = (6.20D-11)
  RCONST(717) = (SUN*5.50E-5+SUN*3.98E-6)
  RCONST(718) = (7.36D-11)
  RCONST(719) = (SUN*5.79E-6+SUN*1.80E-6+SUN*1.80E-6)
  RCONST(720) = (1.00D+06)
  RCONST(721) = (5.22D-11)
  RCONST(722) = (SUN*2.68E-4+SUN*1.80E-6+SUN*1.80E-6)
  RCONST(723) = (9.76D-11)
  RCONST(724) = (SUN*1.80E-6+SUN*1.80E-6)
  RCONST(725) = (GCARR(3.15E-14,0.0E+00,920.0))
  RCONST(726) = (3.13D-13)
  RCONST(727) = (1.53D-12)
  RCONST(728) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(729) = (2.7D-12*EXP(360/TEMP))
! RCONST(730) = constant rate coefficient
  RCONST(731) = (2.50D-13*RO2)
  RCONST(732) = (1.13D-10)
  RCONST(733) = (SUN*5.79E-6)
  RCONST(734) = (1.00D+06*0.85)
  RCONST(735) = (1.00D+06*0.15)
  RCONST(736) = (1.10D-10)
  RCONST(737) = (1.07D-10)
  RCONST(738) = (SUN*2.99E-6)
  RCONST(739) = (SUN*5.79E-6)
  RCONST(740) = (1.00D+06)
  RCONST(741) = (1.03D-10)
  RCONST(742) = (SUN*2.99E-6)
  RCONST(743) = (2.08D-12)
  RCONST(744) = (2.86D-13)
  RCONST(745) = (1.00D+06)
  RCONST(746) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(747) = (5.2D-13*EXP(980/TEMP)*0.15)
  RCONST(748) = (5.2D-13*EXP(980/TEMP)*0.41)
  RCONST(749) = (7.5D-12*EXP(290/TEMP))
  RCONST(750) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(751) = (2.3E-12*1.74)
  RCONST(752) = (1.00D-11*0.70*RO2)
  RCONST(753) = (1.00D-11*0.30*RO2)
  RCONST(754) = (2.91D-13*EXP(1300/TEMP)*0.820)
  RCONST(755) = (2.7D-12*EXP(360/TEMP))
! RCONST(756) = constant rate coefficient
  RCONST(757) = (2.50D-13*RO2)
  RCONST(758) = (5.00D-18)
  RCONST(759) = (4.90D-11)
  RCONST(760) = (SUN*8.92E-3*0.0057)
  RCONST(761) = (5.87D-11)
  RCONST(762) = (SUN*1.80E-6+SUN*1.80E-6)
  RCONST(763) = (6.18D-11)
  RCONST(764) = (SUN*5.79E-6+SUN*1.80E-6+SUN*1.80E-6)
  RCONST(765) = (5.85D-11)
  RCONST(766) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(767) = (7.26D-11)
  RCONST(768) = (SUN*5.79E-6+SUN*5.50E-5+SUN*3.98E-6)
  RCONST(769) = (2.92D-11)
  RCONST(770) = (SUN*5.50E-5+SUN*3.98E-6)
  RCONST(771) = (1.00D+06*0.5)
  RCONST(772) = (1.00D+06*0.5)
  RCONST(773) = (3.58D-11)
  RCONST(774) = (SUN*5.50E-5+SUN*2.68E-4)
  RCONST(775) = (6.03D-11)
  RCONST(776) = (SUN*5.50E-5+SUN*3.98E-6)
  RCONST(777) = (8.25D-11)
  RCONST(778) = (SUN*1.80E-6+SUN*1.80E-6)
  RCONST(779) = (SUN*5.79E-6)
  RCONST(780) = (SUN*5.50E-5)
  RCONST(781) = (1.00D+06)
  RCONST(782) = (1.14D-10)
  RCONST(783) = (SUN*1.80E-6+SUN*1.80E-6)
  RCONST(784) = (SUN*5.50E-5)
  RCONST(785) = (1.43D-10)
  RCONST(786) = (SUN*5.79E-6)
  RCONST(787) = (SUN*2.99E-6)
  RCONST(788) = (1.00D+06)
  RCONST(789) = (1.40D-10)
  RCONST(790) = (SUN*2.99E-6)
  RCONST(791) = (3.13D-13)
  RCONST(792) = (1.53D-12)
  RCONST(793) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(794) = (2.7D-12*EXP(360/TEMP))
! RCONST(795) = constant rate coefficient
  RCONST(796) = (2.50D-13*RO2)
  RCONST(797) = (1.51D-10)
  RCONST(798) = (SUN*5.79E-6)
  RCONST(799) = (1.00D+06*0.68)
  RCONST(800) = (1.00D+06*0.32)
  RCONST(801) = (1.48D-10)
  RCONST(802) = (2.08D-12)
  RCONST(803) = (2.86D-13)
  RCONST(804) = (1.00D+06)
  RCONST(805) = (5.2D-13*EXP(980/TEMP)*0.15)
  RCONST(806) = (5.2D-13*EXP(980/TEMP)*0.41)
  RCONST(807) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(808) = (7.5D-12*EXP(290/TEMP))
  RCONST(809) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(810) = (2.3E-12*1.74)
  RCONST(811) = (1.00D-11*0.30*RO2)
  RCONST(812) = (1.00D-11*0.70*RO2)
  RCONST(813) = (4.06D-11)
  RCONST(814) = (SUN*1.10E-6+SUN*1.10E-6)
  RCONST(815) = (4.37D-11)
  RCONST(816) = (SUN*5.79E-6+SUN*1.10E-6+SUN*1.10E-6)
  RCONST(817) = (4.05D-11)
  RCONST(818) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(819) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(820) = (5.2D-13*EXP(980/TEMP)*0.15)
  RCONST(821) = (5.2D-13*EXP(980/TEMP)*0.41)
  RCONST(822) = (7.5D-12*EXP(290/TEMP))
  RCONST(823) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(824) = (2.3E-12*1.74)
  RCONST(825) = (1.00D-11*RO2*0.7)
  RCONST(826) = (1.00D-11*RO2*0.3)
  RCONST(827) = (2.91D-13*EXP(1300/TEMP)*0.520)
  RCONST(828) = (2.7D-12*EXP(360/TEMP))
! RCONST(829) = constant rate coefficient
  RCONST(830) = (8.80D-13*RO2*0.2)
  RCONST(831) = (8.80D-13*RO2*0.2)
  RCONST(832) = (8.80D-13*RO2*0.6)
  RCONST(833) = (1.31D-10)
  RCONST(834) = (SUN*5.79E-6+SUN*1.99E-5*2)
  RCONST(835) = (4.38D-11)
  RCONST(836) = (SUN*5.50E-5)
  RCONST(837) = (1.00D+06*0.5)
  RCONST(838) = (1.00D+06*0.5)
  RCONST(839) = (9.20D-11)
  RCONST(840) = (SUN*1.99E-5*2+SUN*3.98E-6)
  RCONST(841) = (8.23D-11)
  RCONST(842) = (SUN*1.99E-5*2)
  RCONST(843) = (7.52D-11)
  RCONST(844) = (SUN*5.79E-6)
  RCONST(845) = (SUN*1.10E-6+SUN*1.10E-6)
  RCONST(846) = (1.00D+06)
  RCONST(847) = (4.90D-11)
  RCONST(848) = (SUN*1.25E-4+SUN*1.10E-6+SUN*1.10E-6)
  RCONST(849) = (7.75D-11)
  RCONST(850) = (SUN*1.10E-6+SUN*1.10E-6)
  RCONST(851) = (2.91D-13*EXP(1300/TEMP)*0.625)
  RCONST(852) = (2.7D-12*EXP(360/TEMP))
! RCONST(853) = constant rate coefficient
  RCONST(854) = (8.80D-13*0.20*RO2)
  RCONST(855) = (8.80D-13*0.20*RO2)
  RCONST(856) = (8.80D-13*0.60*RO2)
  RCONST(857) = (9.00D-14)
  RCONST(858) = (9.00D-13)
  RCONST(859) = (2.91D-13*EXP(1300/TEMP)*0.770)
  RCONST(860) = (2.7D-12*EXP(360/TEMP))
! RCONST(861) = constant rate coefficient
  RCONST(862) = (2.50D-13*RO2)
  RCONST(863) = (1.07D-10)
  RCONST(864) = (SUN*2.99E-6)
  RCONST(865) = (SUN*5.79E-6)
  RCONST(866) = (1.00D+06)
  RCONST(867) = (1.04D-10)
  RCONST(868) = (SUN*2.99E-6)
  RCONST(869) = (2.08D-12)
  RCONST(870) = (2.86D-13)
  RCONST(871) = (1.00D+06)
  RCONST(872) = (1.16D-10)
  RCONST(873) = (SUN*5.79E-6)
  RCONST(874) = (1.00D+06*0.71)
  RCONST(875) = (1.00D+06*0.29)
  RCONST(876) = (1.13D-10)
  RCONST(877) = (1.1D-11*EXP(364/TEMP)*0.65)
  RCONST(878) = (1.1D-11*EXP(364/TEMP)*0.20)
  RCONST(879) = (1.1D-11*EXP(364/TEMP)*0.15)
  RCONST(880) = (7.5D-12*EXP(290/TEMP))
  RCONST(881) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(882) = (2.3E-12*1.74)
  RCONST(883) = (1.00D-11*RO2*0.7)
  RCONST(884) = (1.00D-11*RO2*0.3)
  RCONST(885) = (1.15D-10)
  RCONST(886) = (SUN*5.79E-6)
  RCONST(887) = (1.00D+06*0.68)
  RCONST(888) = (1.00D+06*0.32)
  RCONST(889) = (1.11D-10)
  RCONST(890) = (1.07D-10)
  RCONST(891) = (SUN*2.99E-6)
  RCONST(892) = (SUN*5.79E-6)
  RCONST(893) = (1.00D+06)
  RCONST(894) = (1.04D-10)
  RCONST(895) = (SUN*2.99E-6)
  RCONST(896) = (3.13D-13*1.0)
  RCONST(897) = (2.8D-12)
  RCONST(898) = (2.08D-12)
  RCONST(899) = (2.86D-13)
  RCONST(900) = (1.00D+06*1.0)
  RCONST(901) = (5.98D-11)
  RCONST(902) = (SUN*1.10E-6+SUN*1.10E-6)
  RCONST(903) = (6.29D-11)
  RCONST(904) = (SUN*5.79E-6+SUN*1.10E-6+SUN*1.10E-6)
  RCONST(905) = (5.96D-11)
  RCONST(906) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(907) = (7.04D-11)
  RCONST(908) = (SUN*5.79E-6+SUN*5.50E-5+SUN*3.98E-6)
  RCONST(909) = (3.06D-11)
  RCONST(910) = (SUN*5.50E-5+SUN*3.98E-6)
  RCONST(911) = (1.00D+06*0.50)
  RCONST(912) = (1.00D+06*0.50)
  RCONST(913) = (4.06D-11)
  RCONST(914) = (SUN*5.50E-5+SUN*2.68E-4)
  RCONST(915) = (6.17D-11)
  RCONST(916) = (SUN*5.50E-5+SUN*3.98E-6)
  RCONST(917) = (9.42D-11)
  RCONST(918) = (SUN*1.80E-6+SUN*1.80E-6)
  RCONST(919) = (SUN*5.79E-6)
  RCONST(920) = (1.00D+06)
  RCONST(921) = (6.45D-11)
  RCONST(922) = (SUN*1.25E-4+SUN*1.80E-6+SUN*1.80E-6)
  RCONST(923) = (1.18D-10)
  RCONST(924) = (SUN*1.80E-6+SUN*1.80E-6)
  RCONST(925) = (1.10D-12)
  RCONST(926) = (4.66D-12)
  RCONST(927) = (SUN*5.79E-6)
  RCONST(928) = (1.06D-12)
  RCONST(929) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3)*0.67)
  RCONST(930) = (4.65D-11)
  RCONST(931) = (SUN*5.79E-6)
  RCONST(932) = (2.91D-13*EXP(1300/TEMP)*0.625)
  RCONST(933) = (2.7D-12*EXP(360/TEMP))
! RCONST(934) = constant rate coefficient
  RCONST(935) = (8.80D-13*RO2)
  RCONST(936) = (1.00D+06*0.5)
  RCONST(937) = (1.00D+06*0.5)
  RCONST(938) = (2.91D-13*EXP(1300/TEMP)*0.706)
  RCONST(939) = (2.7D-12*EXP(360/TEMP))
! RCONST(940) = constant rate coefficient
  RCONST(941) = (8.80D-13*RO2*0.2)
  RCONST(942) = (8.80D-13*RO2*0.6)
  RCONST(943) = (8.80D-13*RO2*0.2)
  RCONST(944) = (2.91D-13*EXP(1300/TEMP)*0.625)
  RCONST(945) = (2.7D-12*EXP(360/TEMP))
! RCONST(946) = constant rate coefficient
  RCONST(947) = (2.00D-12*0.6*RO2)
  RCONST(948) = (2.00D-12*0.2*RO2)
  RCONST(949) = (2.00D-12*0.2*RO2)
  RCONST(950) = (2.91D-13*EXP(1300/TEMP)*0.706)
  RCONST(951) = (2.7D-12*EXP(360/TEMP))
! RCONST(952) = constant rate coefficient
  RCONST(953) = (1.40D-12*0.20*RO2)
  RCONST(954) = (1.40D-12*0.20*RO2)
  RCONST(955) = (1.40D-12*0.60*RO2)
  RCONST(956) = (5.2D-13*EXP(980/TEMP)*0.15)
  RCONST(957) = (5.2D-13*EXP(980/TEMP)*0.41)
  RCONST(958) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(959) = (7.5D-12*EXP(290/TEMP))
  RCONST(960) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(961) = (2.3E-12*1.74)
  RCONST(962) = (1.00D-11*0.30*RO2)
  RCONST(963) = (1.00D-11*0.70*RO2)
  RCONST(964) = (1.00D+06*0.11)
  RCONST(965) = (1.00D+06*0.89)
  RCONST(966) = (1.00D+06*0.27)
  RCONST(967) = (1.00D+06*0.255)
  RCONST(968) = (1.00D+06*0.255)
  RCONST(969) = (1.00D+06*0.22)
  RCONST(970) = (2.91D-13*EXP(1300/TEMP)*0.770)
  RCONST(971) = (2.7D-12*EXP(360/TEMP))
! RCONST(972) = constant rate coefficient
  RCONST(973) = (9.20D-14*0.3*RO2)
  RCONST(974) = (9.20D-14*0.7*RO2)
  RCONST(975) = (1.00D-12)
  RCONST(976) = (8.00D-17)
  RCONST(977) = (2.66D-11)
  RCONST(978) = (7.00D11*EXP(-3160/TEMP)+5.00D-12)
  RCONST(979) = (5.00D-12*3.2*(1-EXP(-550/TEMP)))
  RCONST(980) = (5.00D-12*3.2*EXP(-550/TEMP))
  RCONST(981) = (2.91D-13*EXP(1300/TEMP)*0.770)
  RCONST(982) = (2.7D-12*EXP(360/TEMP))
! RCONST(983) = constant rate coefficient
  RCONST(984) = (9.20D-14*RO2)
  RCONST(985) = (2.91D-13*EXP(1300/TEMP)*0.770)
  RCONST(986) = (2.7D-12*EXP(360/TEMP))
  RCONST(987) = (2.3E-12*1.74)
  RCONST(988) = (8.80D-13*RO2*0.2)
  RCONST(989) = (8.80D-13*RO2*0.6)
  RCONST(990) = (8.80D-13*RO2*0.2)
  RCONST(991) = (2.91D-13*EXP(1300/TEMP)*0.770)
  RCONST(992) = (2.7D-12*EXP(360/TEMP))
! RCONST(993) = constant rate coefficient
  RCONST(994) = (9.20D-14*RO2*0.70)
  RCONST(995) = (9.20D-14*RO2*0.30)
  RCONST(996) = (2.91D-13*EXP(1300/TEMP)*0.706)
  RCONST(997) = (2.7D-12*EXP(360/TEMP))
! RCONST(998) = constant rate coefficient
  RCONST(999) = (9.20D-14*RO2)
  RCONST(1000) = (1.00D+06*0.5)
  RCONST(1001) = (1.00D+06*0.5)
  RCONST(1002) = (2.91D-13*EXP(1300/TEMP)*0.706)
  RCONST(1003) = (2.7D-12*EXP(360/TEMP))
! RCONST(1004) = constant rate coefficient
  RCONST(1005) = (9.20D-14*RO2*0.7)
  RCONST(1006) = (9.20D-14*RO2*0.3)
  RCONST(1007) = (2.22D-11)
  RCONST(1008) = (SUN*5.50E-5)
  RCONST(1009) = (2.53D-11)
  RCONST(1010) = (SUN*5.79E-6+SUN*5.50E-5)
  RCONST(1011) = (2.21D-11)
  RCONST(1012) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(1013) = (6.97D-11)
  RCONST(1014) = (SUN*5.50E-5)
  RCONST(1015) = (SUN*5.79E-6)
  RCONST(1016) = (SUN*3.98E-6)
  RCONST(1017) = (1.00D+06*0.5)
  RCONST(1018) = (1.00D+06*0.5)
  RCONST(1019) = (SUN*1.25E-4)
  RCONST(1020) = (1.23D-11)
  RCONST(1021) = (1.36D-10)
  RCONST(1022) = (SUN*1.99E-5*2)
  RCONST(1023) = (6.63D-11)
  RCONST(1024) = (SUN*5.50E-5)
  RCONST(1025) = (SUN*3.98E-6)
  RCONST(1026) = (5.08D-11)
  RCONST(1027) = (SUN*1.80E-6)
  RCONST(1028) = (2.08D-12)
  RCONST(1029) = (2.86D-13)
  RCONST(1030) = (4.65D-11)
  RCONST(1031) = (SUN*5.79E-6)
  RCONST(1032) = (1.00D-12)
  RCONST(1033) = (2.35D-10)
  RCONST(1034) = (5.03D-12)
  RCONST(1035) = (6.83D-12)
  RCONST(1036) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(1037) = (2.7D-12*EXP(360/TEMP))
! RCONST(1038) = constant rate coefficient
  RCONST(1039) = (8.80D-13*RO2)
  RCONST(1040) = (5.44D-11)
  RCONST(1041) = (SUN*1.80E-6+SUN*1.80E-6)
  RCONST(1042) = (SUN*1.25E-4)
  RCONST(1043) = (8.0D-13)
  RCONST(1044) = (1.10D-12)
  RCONST(1045) = (4.66D-12)
  RCONST(1046) = (SUN*5.79E-6)
  RCONST(1047) = (1.06D-12)
  RCONST(1048) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3)*0.67)
  RCONST(1049) = (1.90D-12*EXP(190/TEMP))
  RCONST(1050) = (SUN*5.79E-6)
  RCONST(1051) = (3.90D-13)
  RCONST(1052) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(1053) = (5.2D-13*EXP(980/TEMP)*0.56)
  RCONST(1054) = (7.5D-12*EXP(290/TEMP))
  RCONST(1055) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(1056) = (2.3E-12*1.74)
  RCONST(1057) = (1.00D-11*0.70*RO2)
  RCONST(1058) = (1.00D-11*0.30*RO2)
  RCONST(1059) = (1.00D+06*0.89)
  RCONST(1060) = (1.00D+06*0.11)
  RCONST(1061) = (1.00D+06*0.11)
  RCONST(1062) = (1.00D+06*0.89)
  RCONST(1063) = (2.91D-13*EXP(1300/TEMP)*0.770)
  RCONST(1064) = (2.7D-12*EXP(360/TEMP))
! RCONST(1065) = constant rate coefficient
  RCONST(1066) = (9.20D-14*0.70*RO2)
  RCONST(1067) = (9.20D-14*0.30*RO2)
  RCONST(1068) = (2.91D-13*EXP(1300/TEMP)*0.770)
  RCONST(1069) = (2.7D-12*EXP(360/TEMP))
! RCONST(1070) = constant rate coefficient
  RCONST(1071) = (9.20D-14*RO2)
  RCONST(1072) = (2.91D-13*EXP(1300/TEMP)*0.706)
  RCONST(1073) = (2.7D-12*EXP(360/TEMP))
! RCONST(1074) = constant rate coefficient
  RCONST(1075) = (8.80D-13*RO2*0.2)
  RCONST(1076) = (8.80D-13*RO2*0.6)
  RCONST(1077) = (8.80D-13*RO2*0.2)
  RCONST(1078) = (2.91D-13*EXP(1300/TEMP)*0.770)
  RCONST(1079) = (2.7D-12*EXP(360/TEMP))
! RCONST(1080) = constant rate coefficient
  RCONST(1081) = (9.20D-14*RO2*0.70)
  RCONST(1082) = (9.20D-14*RO2*0.30)
  RCONST(1083) = (5.2D-13*EXP(980/TEMP)*0.41)
  RCONST(1084) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(1085) = (5.2D-13*EXP(980/TEMP)*0.15)
  RCONST(1086) = (7.5D-12*EXP(290/TEMP))
  RCONST(1087) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(1088) = (2.3E-12*1.74)
  RCONST(1089) = (1.00D-11*0.70*RO2)
  RCONST(1090) = (1.00D-11*0.30*RO2)
  RCONST(1091) = (2.91D-13*EXP(1300/TEMP)*0.706)
  RCONST(1092) = (2.7D-12*EXP(360/TEMP))
! RCONST(1093) = constant rate coefficient
  RCONST(1094) = (9.20D-14*RO2*0.30)
  RCONST(1095) = (9.20D-14*RO2*0.70)
  RCONST(1096) = (2.91D-13*EXP(1300/TEMP)*0.706)
  RCONST(1097) = (2.7D-12*EXP(360/TEMP))
! RCONST(1098) = constant rate coefficient
  RCONST(1099) = (9.20D-14*RO2)
  RCONST(1100) = (2.91D-13*EXP(1300/TEMP)*0.625)
  RCONST(1101) = (2.7D-12*EXP(360/TEMP))
! RCONST(1102) = constant rate coefficient
  RCONST(1103) = (2.00D-12*RO2)
  RCONST(1104) = (2.91D-13*EXP(1300/TEMP)*0.706)
  RCONST(1105) = (2.7D-12*EXP(360/TEMP))
! RCONST(1106) = constant rate coefficient
  RCONST(1107) = (9.20D-14*RO2*0.70)
  RCONST(1108) = (9.20D-14*RO2*0.30)
  RCONST(1109) = (1.00D+06*0.35)
  RCONST(1110) = (1.00D+06*0.65)
  RCONST(1111) = (2.40D-17*0.27)
  RCONST(1112) = (2.40D-17*0.73)
  RCONST(1113) = (4.72D-11)
  RCONST(1114) = (SUN*5.64E-4*2)
  RCONST(1115) = (4.37D-11)
  RCONST(1116) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(1117) = (1.00D+06*0.83)
  RCONST(1118) = (1.00D+06*0.17)
  RCONST(1119) = (3.59D-12)
  RCONST(1120) = (SUN*5.64E-4*2)
  RCONST(1121) = (5.43D-11)
  RCONST(1122) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(1123) = (2.27D-11)
  RCONST(1124) = (SUN*5.50E-5)
  RCONST(1125) = (2.57D-11)
  RCONST(1126) = (SUN*5.79E-6+SUN*5.50E-5)
  RCONST(1127) = (2.25D-11)
  RCONST(1128) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(1129) = (1.35D-10)
  RCONST(1130) = (SUN*5.79E-6)
  RCONST(1131) = (SUN*5.50E-5*2)
  RCONST(1132) = (1.00D+06)
  RCONST(1133) = (4.18D-11)
  RCONST(1134) = (SUN*5.50E-5*2)
  RCONST(1135) = (2.91D-13*EXP(1300/TEMP)*0.387)
  RCONST(1136) = (2.7D-12*EXP(360/TEMP))
! RCONST(1137) = constant rate coefficient
  RCONST(1138) = (2.00D-12*0.2*RO2)
  RCONST(1139) = (2.00D-12*0.6*RO2)
  RCONST(1140) = (2.00D-12*0.2*RO2)
  RCONST(1141) = (SUN*1.99E-5)
  RCONST(1142) = (1.4D-12*EXP(-1860/TEMP)*4.0)
  RCONST(1143) = (2.45D-11)
  RCONST(1144) = (2.08D-12)
  RCONST(1145) = (2.86D-13)
  RCONST(1146) = (4.65D-11)
  RCONST(1147) = (SUN*5.79E-6)
  RCONST(1148) = (1.00D-12)
  RCONST(1149) = (2.35D-11)
  RCONST(1150) = (5.03D-12)
  RCONST(1151) = (6.83D-12)
  RCONST(1152) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(1153) = (2.7D-12*EXP(360/TEMP))
! RCONST(1154) = constant rate coefficient
  RCONST(1155) = (8.80D-13*RO2)
  RCONST(1156) = (4.66D-11)
  RCONST(1157) = (SUN*1.10E-6+SUN*1.10E-6)
  RCONST(1158) = (1.10D-12)
  RCONST(1159) = (4.66D-12)
  RCONST(1160) = (SUN*5.79E-6)
  RCONST(1161) = (1.06D-12)
  RCONST(1162) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3)*0.67)
  RCONST(1163) = (1.90D-12*EXP(190/TEMP))
  RCONST(1164) = (SUN*5.79E-6)
  RCONST(1165) = (3.90D-13)
  RCONST(1166) = (2.91D-13*EXP(1300/TEMP)*0.770)
  RCONST(1167) = (2.7D-12*EXP(360/TEMP))
! RCONST(1168) = constant rate coefficient
  RCONST(1169) = (8.80D-13*0.20*RO2)
  RCONST(1170) = (8.80D-13*0.20*RO2)
  RCONST(1171) = (8.80D-13*0.60*RO2)
  RCONST(1172) = (2.91D-13*EXP(1300/TEMP)*0.625)
  RCONST(1173) = (2.7D-12*EXP(360/TEMP))
! RCONST(1174) = constant rate coefficient
  RCONST(1175) = (8.80D-12*RO2)
  RCONST(1176) = (2.08D-12)
  RCONST(1177) = (2.86D-13)
  RCONST(1178) = (4.65D-11)
  RCONST(1179) = (SUN*5.79E-6)
  RCONST(1180) = (1.00D-12)
  RCONST(1181) = (2.35D-11)
  RCONST(1182) = (5.03D-12)
  RCONST(1183) = (6.83D-12)
  RCONST(1184) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(1185) = (2.7D-12*EXP(360/TEMP))
! RCONST(1186) = constant rate coefficient
  RCONST(1187) = (8.80D-13*RO2)
  RCONST(1188) = (4.38D-11)
  RCONST(1189) = (SUN*1.10E-6+SUN*1.10E-6)
  RCONST(1190) = (3.70D-11)
  RCONST(1191) = (SUN*1.10E-6+SUN*1.10E-6)
  RCONST(1192) = (4.00D-11)
  RCONST(1193) = (SUN*5.64E-4*2)
  RCONST(1194) = (1.00D+06*0.40)
  RCONST(1195) = (1.00D+06*0.60)
  RCONST(1196) = (3.70D-11)
  RCONST(1197) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(1198) = (2.31D-11)
  RCONST(1199) = (SUN*5.50E-5)
  RCONST(1200) = (2.62D-11)
  RCONST(1201) = (SUN*5.79E-6+SUN*5.50E-5)
  RCONST(1202) = (2.29D-11)
  RCONST(1203) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(1204) = (1.44D-10)
  RCONST(1205) = (SUN*5.79E-6)
  RCONST(1206) = (SUN*1.99E-5*2)
  RCONST(1207) = (1.00D+06)
  RCONST(1208) = (5.77D-11)
  RCONST(1209) = (SUN*1.99E-5*2)
  RCONST(1210) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(1211) = (5.2D-13*EXP(980/TEMP)*0.56)
  RCONST(1212) = (7.5D-12*EXP(290/TEMP))
  RCONST(1213) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(1214) = (2.3E-12*1.74)
  RCONST(1215) = (1.00D-11*RO2)
  RCONST(1216) = (1.22D-10)
  RCONST(1217) = (1.90D-12*EXP(190/TEMP))
  RCONST(1218) = (SUN*1.99E-5*2)
  RCONST(1219) = (SUN*5.79E-6)
  RCONST(1220) = (1.00D+06)
  RCONST(1221) = (3.67D-11)
  RCONST(1222) = (SUN*1.25E-4+SUN*1.99E-5)
  RCONST(1223) = (8.13D-11)
  RCONST(1224) = (SUN*1.99E-5*2)
  RCONST(1225) = (2.08D-12)
  RCONST(1226) = (2.86D-13)
  RCONST(1227) = (3.60D-12)
  RCONST(1228) = (SUN*5.79E-6)
  RCONST(1229) = (2.60D-12)
  RCONST(1230) = (3.47D-12)
  RCONST(1231) = (2.91D-13*EXP(1300/TEMP)*0.770)
  RCONST(1232) = (2.7D-12*EXP(360/TEMP))
! RCONST(1233) = constant rate coefficient
  RCONST(1234) = (8.80D-13*RO2)
  RCONST(1235) = (3.00D-13)
  RCONST(1236) = (4.6D-12)
  RCONST(1237) = (4.66D-12)
  RCONST(1238) = (SUN*5.79E-6)
  RCONST(1239) = (1.10D-12)
  RCONST(1240) = (1.06D-12)
  RCONST(1241) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3)*0.67)
  RCONST(1242) = (1.00D-12)
  RCONST(1243) = (2.3D-11)
  RCONST(1244) = (2.08D-12)
  RCONST(1245) = (2.86D-13)
  RCONST(1246) = (5.03D-12)
  RCONST(1247) = (6.83D-12)
  RCONST(1248) = (2.91D-13*EXP(1300/TEMP)*0.820)
  RCONST(1249) = (2.7D-12*EXP(360/TEMP))
! RCONST(1250) = constant rate coefficient
  RCONST(1251) = (8.80D-13*RO2)
  RCONST(1252) = (6.18D-12)
  RCONST(1253) = (SUN*5.79E-6)
  RCONST(1254) = (1.00D+06*0.5)
  RCONST(1255) = (1.00D+06*0.5)
  RCONST(1256) = (1.20D-15)
  RCONST(1257) = (1.00D-14)
  RCONST(1258) = (1.00D-15)
  RCONST(1259) = (7.00D-14)
  RCONST(1260) = (6.00D-18)
  RCONST(1261) = (1.00D-17)
  RCONST(1262) = (3.68D-11)
  RCONST(1263) = (SUN*5.79E-6)
  RCONST(1264) = (1.00D+06)
  RCONST(1265) = (1.78D-11)
  RCONST(1266) = (2.78D-11)
  RCONST(1267) = (SUN*5.79E-6+SUN*2.68E-4)
  RCONST(1268) = (1.90D-12*EXP(190/TEMP))
  RCONST(1269) = (5.99D-12)
  RCONST(1270) = (1.00D+06)
  RCONST(1271) = (SUN*2.68E-4)
  RCONST(1272) = (2.69D-12)
  RCONST(1273) = (SUN*1.25E-4)
  RCONST(1274) = (SUN*2.68E-4)
  RCONST(1275) = (1.4D-12*EXP(-1860/TEMP)*4.0)
  RCONST(1276) = (1.23D-11)
  RCONST(1277) = (1.21D-10)
  RCONST(1278) = (1.90D-12*EXP(190/TEMP))
  RCONST(1279) = (SUN*5.79E-6+SUN*1.99E-5+SUN*3.98E-6)
  RCONST(1280) = (1.00D+06)
  RCONST(1281) = (7.48D-11)
  RCONST(1282) = (SUN*1.99E-5+SUN*2.68E-4)
  RCONST(1283) = (6.26D-11)
  RCONST(1284) = (SUN*1.99E-5+SUN*3.98E-6)
  RCONST(1285) = (5.44D-11)
  RCONST(1286) = (SUN*5.64E-4*2)
  RCONST(1287) = (1.00D+06*0.35)
  RCONST(1288) = (1.00D+06*0.65)
  RCONST(1289) = (5.08D-11)
  RCONST(1290) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(1291) = (8.8D-12*EXP(-1320/TEMP)+1.7D-14*EXP(423/TEMP))
  RCONST(1292) = (SUN*6.09E-7)
  RCONST(1293) = (1.20D-15)
  RCONST(1294) = (1.00D-14)
  RCONST(1295) = (1.00D-15)
  RCONST(1296) = (7.00D-14)
  RCONST(1297) = (6.00D-18)
  RCONST(1298) = (3.12D-11)
  RCONST(1299) = (SUN*3.98E-6+SUN*5.50E-5)
  RCONST(1300) = (SUN*5.79E-6)
  RCONST(1301) = (1.00D+06)
  RCONST(1302) = (2.77D-11)
  RCONST(1303) = (SUN*3.98E-6+SUN*5.50E-5)
  RCONST(1304) = (2.91D-13*EXP(1300/TEMP)*0.706)
  RCONST(1305) = (2.7D-12*EXP(360/TEMP))
! RCONST(1306) = constant rate coefficient
  RCONST(1307) = (9.20D-14*RO2)
  RCONST(1308) = (1.85D-11)
  RCONST(1309) = (SUN*3.98E-6+SUN*1.25E-4)
  RCONST(1310) = (2.91D-13*EXP(1300/TEMP)*0.770)
  RCONST(1311) = (2.7D-12*EXP(360/TEMP))
! RCONST(1312) = constant rate coefficient
  RCONST(1313) = (9.20D-14*RO2*0.70)
  RCONST(1314) = (9.20D-14*RO2*0.30)
  RCONST(1315) = (5.2D-13*EXP(980/TEMP)*0.15)
  RCONST(1316) = (5.2D-13*EXP(980/TEMP)*0.41)
  RCONST(1317) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(1318) = (7.5D-12*EXP(290/TEMP))
  RCONST(1319) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(1320) = (2.3E-12*1.74)
  RCONST(1321) = (1.00D-11*0.7*RO2)
  RCONST(1322) = (1.00D-11*0.3*RO2)
  RCONST(1323) = (4.95D-12)
  RCONST(1324) = (SUN*5.79E-6)
  RCONST(1325) = (1.00D+06)
  RCONST(1326) = (4.42D-11)
  RCONST(1327) = (1.90D-12*EXP(190/TEMP))
  RCONST(1328) = (SUN*5.79E-6)
  RCONST(1329) = (1.00D+06)
  RCONST(1330) = (4.65D-13)
  RCONST(1331) = (SUN*3.98E-6)
  RCONST(1332) = (1.91D-11)
  RCONST(1333) = (7.12D-12)
  RCONST(1334) = (SUN*5.79E-6)
  RCONST(1335) = (1.00D+06)
  RCONST(1336) = (3.76D-12)
  RCONST(1337) = (4.44D-12)
  RCONST(1338) = (SUN*5.79E-6)
  RCONST(1339) = (1.00D+06)
  RCONST(1340) = (2.91D-13*EXP(1300/TEMP)*0.625)
  RCONST(1341) = (2.7D-12*EXP(360/TEMP))
! RCONST(1342) = constant rate coefficient
  RCONST(1343) = (2.00D-12*RO2)
  RCONST(1344) = (1.20D-15)
  RCONST(1345) = (1.00D-14)
  RCONST(1346) = (1.00D-15)
  RCONST(1347) = (7.00D-14)
  RCONST(1348) = (6.00D-18)
  RCONST(1349) = (1.00D-17)
  RCONST(1350) = (2.53D-11)
  RCONST(1351) = (SUN*5.79E-6)
  RCONST(1352) = (1.00D+06)
  RCONST(1353) = (2.19D-11)
  RCONST(1354) = (7.83D-15)
  RCONST(1355) = (5.10D-14)
  RCONST(1356) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(1357) = (2.7D-12*EXP(360/TEMP))
! RCONST(1358) = constant rate coefficient
  RCONST(1359) = (2.50D-13*RO2)
  RCONST(1360) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(1361) = (2.7D-12*EXP(360/TEMP))
! RCONST(1362) = constant rate coefficient
  RCONST(1363) = (8.80D-13*RO2)
  RCONST(1364) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(1365) = (2.7D-12*EXP(360/TEMP))
! RCONST(1366) = constant rate coefficient
  RCONST(1367) = (8.80D-13*RO2*0.2)
  RCONST(1368) = (8.80D-13*RO2*0.6)
  RCONST(1369) = (8.80D-13*RO2*0.2)
  RCONST(1370) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(1371) = (2.7D-12*EXP(360/TEMP))
! RCONST(1372) = constant rate coefficient
  RCONST(1373) = (8.00D-13*RO2)
  RCONST(1374) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(1375) = (2.7D-12*EXP(360/TEMP))
! RCONST(1376) = constant rate coefficient
  RCONST(1377) = (8.00D-13*RO2)
  RCONST(1378) = (2.05D-10)
  RCONST(1379) = (SUN*5.79E-6)
  RCONST(1380) = (1.00D+06*0.35)
  RCONST(1381) = (1.00D+06*0.65)
  RCONST(1382) = (3.90D-11)
  RCONST(1383) = (SUN*5.64E-4*2)
  RCONST(1384) = (3.54D-11)
  RCONST(1385) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(1386) = (3.54D-11)
  RCONST(1387) = (SUN*1.80E-6)
  RCONST(1388) = (1.2D-15)
  RCONST(1389) = (1.0D-14)
  RCONST(1390) = (1.0D-15)
  RCONST(1391) = (7.0D-14)
  RCONST(1392) = (6.0D-18)
  RCONST(1393) = (4.46D-11)
  RCONST(1394) = (SUN*5.50E-5)
  RCONST(1395) = (SUN*5.79E-6+SUN*3.98E-6)
  RCONST(1396) = (1.00D+06)
  RCONST(1397) = (4.13D-11)
  RCONST(1398) = (SUN*3.98E-6+SUN*5.50E-5)
  RCONST(1399) = (5.80D-12)
  RCONST(1400) = (SUN*5.79E-6)
  RCONST(1401) = (1.00D+06)
  RCONST(1402) = (6.01D-11)
  RCONST(1403) = (1.90D-12*EXP(190/TEMP))
  RCONST(1404) = (SUN*5.79E-6)
  RCONST(1405) = (1.00D+06)
  RCONST(1406) = (1.09D-13)
  RCONST(1407) = (SUN*3.98E-6)
  RCONST(1408) = (2.54D-11)
  RCONST(1409) = (3.09D-11)
  RCONST(1410) = (SUN*5.79E-6)
  RCONST(1411) = (1.00D+06)
  RCONST(1412) = (2.74D-11)
  RCONST(1413) = (2.40D-17*0.73)
  RCONST(1414) = (2.40D-17*0.27)
  RCONST(1415) = (4.72D-11)
  RCONST(1416) = (SUN*5.64E-4*2)
  RCONST(1417) = (1.00D+06*0.35)
  RCONST(1418) = (1.00D+06*0.65)
  RCONST(1419) = (4.37D-11)
  RCONST(1420) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(1421) = (4.58D-11)
  RCONST(1422) = (SUN*5.79E-6)
  RCONST(1423) = (SUN*5.50E-5*2)
  RCONST(1424) = (1.00D+06)
  RCONST(1425) = (5.70D-11)
  RCONST(1426) = (SUN*5.50E-5*2)
  RCONST(1427) = (5.16D-12)
  RCONST(1428) = (SUN*5.79E-6)
  RCONST(1429) = (1.00D+06)
  RCONST(1430) = (2.00D-11)
  RCONST(1431) = (SUN*3.98E-6)
  RCONST(1432) = (SUN*5.79E-6)
  RCONST(1433) = (1.00D+06)
  RCONST(1434) = (2.78D-11)
  RCONST(1435) = (SUN*5.79E-6)
  RCONST(1436) = (1.00D+06)
  RCONST(1437) = (2.44D-11)
  RCONST(1438) = (1.50D-12)
  RCONST(1439) = (3.69D-12)
  RCONST(1440) = (SUN*5.79E-6+SUN*2.68E-4)
  RCONST(1441) = (1.00D+06)
  RCONST(1442) = (2.91D-13*EXP(1300/TEMP)*0.706)
  RCONST(1443) = (2.7D-12*EXP(360/TEMP))
! RCONST(1444) = constant rate coefficient
  RCONST(1445) = (9.20D-14*RO2)
  RCONST(1446) = (1.4D-12)
  RCONST(1447) = (5.2D-13*EXP(980/TEMP)*0.56)
  RCONST(1448) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(1449) = (7.5D-12*EXP(290/TEMP))
  RCONST(1450) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(1451) = (2.3E-12*1.74)
  RCONST(1452) = (1.00D-11*RO2)
  RCONST(1453) = (2.91D-11)
  RCONST(1454) = (1.90D-12*EXP(190/TEMP))
  RCONST(1455) = (SUN*5.79E-6)
  RCONST(1456) = (SUN*1.99E-5)
  RCONST(1457) = (1.00D+06)
  RCONST(1458) = (5.2D-13*EXP(980/TEMP)*0.56)
  RCONST(1459) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(1460) = (7.5D-12*EXP(290/TEMP))
  RCONST(1461) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(1462) = (2.3E-12*1.74)
  RCONST(1463) = (1.00D-11*RO2)
  RCONST(1464) = (7.83D-15)
  RCONST(1465) = (5.10D-14)
  RCONST(1466) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(1467) = (2.7D-12*EXP(360/TEMP))
! RCONST(1468) = constant rate coefficient
  RCONST(1469) = (2.50D-13*RO2)
  RCONST(1470) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(1471) = (2.7D-12*EXP(360/TEMP))
! RCONST(1472) = constant rate coefficient
  RCONST(1473) = (9.20D-14*RO2)
  RCONST(1474) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(1475) = (2.7D-12*EXP(360/TEMP))
! RCONST(1476) = constant rate coefficient
  RCONST(1477) = (9.20D-14*RO2*0.7)
  RCONST(1478) = (9.20D-14*RO2*0.3)
  RCONST(1479) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(1480) = (2.7D-12*EXP(360/TEMP))
! RCONST(1481) = constant rate coefficient
  RCONST(1482) = (8.00D-13*RO2)
  RCONST(1483) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(1484) = (2.7D-12*EXP(360/TEMP))
! RCONST(1485) = constant rate coefficient
  RCONST(1486) = (8.00D-13*RO2)
  RCONST(1487) = (2.05D-10)
  RCONST(1488) = (SUN*5.79E-6)
  RCONST(1489) = (3.6D-12)
  RCONST(1490) = (9.75D-11)
  RCONST(1491) = (1.90D-12*EXP(190/TEMP))
  RCONST(1492) = (SUN*3.98E-6*2)
  RCONST(1493) = (SUN*5.79E-6)
  RCONST(1494) = (1.00D+06*0.7)
  RCONST(1495) = (1.00D+06*0.3)
  RCONST(1496) = (4.00D-11)
  RCONST(1497) = (SUN*3.98E-6*2)
  RCONST(1498) = (9.04D-12)
  RCONST(1499) = (SUN*3.98E-6)
  RCONST(1500) = (SUN*2.68E-4)
  RCONST(1501) = (7.83D-11)
  RCONST(1502) = (SUN*5.79E-6)
  RCONST(1503) = (SUN*3.98E-6)
  RCONST(1504) = (SUN*5.50E-5)
  RCONST(1505) = (1.00D+06*0.5)
  RCONST(1506) = (1.00D+06*0.5)
  RCONST(1507) = (7.83D-15)
  RCONST(1508) = (5.10D-14)
  RCONST(1509) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(1510) = (2.7D-12*EXP(360/TEMP))
! RCONST(1511) = constant rate coefficient
  RCONST(1512) = (2.50D-13*RO2)
  RCONST(1513) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(1514) = (2.7D-12*EXP(360/TEMP))
! RCONST(1515) = constant rate coefficient
  RCONST(1516) = (9.20D-14*RO2)
  RCONST(1517) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(1518) = (2.7D-12*EXP(360/TEMP))
! RCONST(1519) = constant rate coefficient
  RCONST(1520) = (9.20D-14*RO2*0.70)
  RCONST(1521) = (9.20D-14*RO2*0.30)
  RCONST(1522) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(1523) = (2.7D-12*EXP(360/TEMP))
! RCONST(1524) = constant rate coefficient
  RCONST(1525) = (8.00D-13*RO2)
  RCONST(1526) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(1527) = (2.7D-12*EXP(360/TEMP))
! RCONST(1528) = constant rate coefficient
  RCONST(1529) = (8.00D-13*RO2)
  RCONST(1530) = (2.05D-10)
  RCONST(1531) = (SUN*5.79E-6)
  RCONST(1532) = (7.33D-11)
  RCONST(1533) = (SUN*5.79E-6)
  RCONST(1534) = (6.97D-11)
  RCONST(1535) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(1536) = (2.45D-11)
  RCONST(1537) = (SUN*1.25E-4*2)
  RCONST(1538) = (2.25D-15)
  RCONST(1539) = (3.00D-14)
  RCONST(1540) = (2.91D-13*EXP(1300/TEMP)*0.770)
  RCONST(1541) = (2.7D-12*EXP(360/TEMP))
! RCONST(1542) = constant rate coefficient
  RCONST(1543) = (2.50D-13*RO2)
  RCONST(1544) = (2.91D-13*EXP(1300/TEMP)*0.770)
  RCONST(1545) = (2.7D-12*EXP(360/TEMP))
! RCONST(1546) = constant rate coefficient
  RCONST(1547) = (8.00D-13*RO2)
  RCONST(1548) = (2.91D-13*EXP(1300/TEMP)*0.770)
  RCONST(1549) = (2.7D-12*EXP(360/TEMP))
! RCONST(1550) = constant rate coefficient
  RCONST(1551) = (8.00D-13*RO2)
  RCONST(1552) = (1.90D-12*EXP(190/TEMP))
  RCONST(1553) = (SUN*5.79E-6)
  RCONST(1554) = (2.91D-13*EXP(1300/TEMP)*0.770)
  RCONST(1555) = (2.7D-12*EXP(360/TEMP))
! RCONST(1556) = constant rate coefficient
  RCONST(1557) = (8.80D-13*RO2)
  RCONST(1558) = (2.91D-13*EXP(1300/TEMP)*0.770)
  RCONST(1559) = (2.7D-12*EXP(360/TEMP))
! RCONST(1560) = constant rate coefficient
  RCONST(1561) = (8.80D-13*RO2*0.2)
  RCONST(1562) = (8.80D-13*RO2*0.6)
  RCONST(1563) = (8.80D-13*RO2*0.2)
  RCONST(1564) = (2.91D-13*EXP(1300/TEMP)*0.820)
  RCONST(1565) = (2.7D-12*EXP(360/TEMP))
! RCONST(1566) = constant rate coefficient
  RCONST(1567) = (8.80D-13*RO2)
  RCONST(1568) = (2.91D-13*EXP(1300/TEMP)*0.820)
  RCONST(1569) = (2.7D-12*EXP(360/TEMP))
! RCONST(1570) = constant rate coefficient
  RCONST(1571) = (8.80D-13*RO2*0.2)
  RCONST(1572) = (8.80D-13*RO2*0.6)
  RCONST(1573) = (8.80D-13*RO2*0.2)
  RCONST(1574) = (7.83D-15)
  RCONST(1575) = (5.10D-14)
  RCONST(1576) = (2.91D-13*EXP(1300/TEMP)*0.820)
  RCONST(1577) = (2.7D-12*EXP(360/TEMP))
! RCONST(1578) = constant rate coefficient
  RCONST(1579) = (2.50D-13*RO2)
  RCONST(1580) = (2.91D-13*EXP(1300/TEMP)*0.820)
  RCONST(1581) = (2.7D-12*EXP(360/TEMP))
! RCONST(1582) = constant rate coefficient
  RCONST(1583) = (8.00D-13*RO2)
  RCONST(1584) = (2.91D-13*EXP(1300/TEMP)*0.820)
  RCONST(1585) = (2.7D-12*EXP(360/TEMP))
! RCONST(1586) = constant rate coefficient
  RCONST(1587) = (8.00D-13*RO2)
  RCONST(1588) = (2.05D-10)
  RCONST(1589) = (SUN*5.79E-6)
  RCONST(1590) = (1.4D-12*EXP(-1860/TEMP)*8.0)
  RCONST(1591) = (3.44D-11)
  RCONST(1592) = (SUN*1.25E-4)
  RCONST(1593) = (1.16D-12)
  RCONST(1594) = (2.19D-11)
  RCONST(1595) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(1596) = (5.2D-13*EXP(980/TEMP)*0.15)
  RCONST(1597) = (5.2D-13*EXP(980/TEMP)*0.41)
  RCONST(1598) = (7.5D-12*EXP(290/TEMP))
  RCONST(1599) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(1600) = (2.3E-12*1.74)
  RCONST(1601) = (1.00D-11*0.7*RO2)
  RCONST(1602) = (1.00D-11*0.3*RO2)
  RCONST(1603) = (1.4D-12*EXP(-1860/TEMP)*5.5)
  RCONST(1604) = (1.72D-11)
  RCONST(1605) = (SUN*1.25E-4+SUN*2.68E-4*2)
  RCONST(1606) = (1.36D-13*EXP(1250/TEMP)*0.15)
  RCONST(1607) = (1.36D-13*EXP(1250/TEMP)*0.85)
  RCONST(1608) = (2.7D-12*EXP(360/TEMP))
! RCONST(1609) = constant rate coefficient
  RCONST(1610) = (2*(3.5D-13*8.0D-12)**0.5*RO2*0.2)
  RCONST(1611) = (2*(3.5D-13*8.0D-12)**0.5*RO2*0.6)
  RCONST(1612) = (2*(3.5D-13*8.0D-12)**0.5*RO2*0.2)
  RCONST(1613) = (9.93D-11)
  RCONST(1614) = (SUN*5.79E-6)
  RCONST(1615) = (1.00D+06)
  RCONST(1616) = (2.91D-13*EXP(1300/TEMP)*0.625)
  RCONST(1617) = (2.7D-12*EXP(360/TEMP))
! RCONST(1618) = constant rate coefficient
  RCONST(1619) = (8.80D-13*0.2*RO2)
  RCONST(1620) = (8.80D-13*0.2*RO2)
  RCONST(1621) = (8.80D-13*0.6*RO2)
  RCONST(1622) = (1.37D-11)
  RCONST(1623) = (SUN*5.79E-6)
  RCONST(1624) = (1.00D+06)
  RCONST(1625) = (1.02D-11)
  RCONST(1626) = (SUN*5.79E-6+SUN*1.99E-5)
  RCONST(1627) = (1.58D-11)
  RCONST(1628) = (2.90D-12)
  RCONST(1629) = (SUN*3.98E-6)
  RCONST(1630) = (2.91D-13*EXP(1300/TEMP))
  RCONST(1631) = (2.7D-12*EXP(360/TEMP))
! RCONST(1632) = constant rate coefficient
  RCONST(1633) = (2.00D-12*RO2)
  RCONST(1634) = (1.4D-12*EXP(-1860/TEMP)*5.5)
  RCONST(1635) = (7.09D-11)
  RCONST(1636) = (SUN*1.99E-5)
  RCONST(1637) = (3.59D-12)
  RCONST(1638) = (SUN*5.79E-6)
  RCONST(1639) = (1.00D+06)
  RCONST(1640) = (8.76D-13)
  RCONST(1641) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(1642) = (2.7D-12*EXP(360/TEMP))
! RCONST(1643) = constant rate coefficient
  RCONST(1644) = (8.00D-13*RO2)
  RCONST(1645) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(1646) = (2.7D-12*EXP(360/TEMP))
! RCONST(1647) = constant rate coefficient
  RCONST(1648) = (8.00D-13*RO2)
  RCONST(1649) = (1.53D-12)
  RCONST(1650) = (SUN*5.79E-6)
  RCONST(1651) = (1.05D-10)
  RCONST(1652) = (SUN*5.79E-6)
  RCONST(1653) = (1.00D+06)
  RCONST(1654) = (1.61D-10)
  RCONST(1655) = (SUN*5.79E-6)
  RCONST(1656) = (1.00D+06)
  RCONST(1657) = (9.85D-11)
  RCONST(1658) = (1.30D-10)
  RCONST(1659) = (1.90D-12*EXP(190/TEMP))
  RCONST(1660) = (SUN*5.79E-6)
  RCONST(1661) = (1.00D+06)
  RCONST(1662) = (1.90D-12*EXP(190/TEMP))
  RCONST(1663) = (SUN*5.79E-6)
  RCONST(1664) = (1.00D+06)
  RCONST(1665) = (4.06D-12)
  RCONST(1666) = (SUN*3.98E-6*2)
  RCONST(1667) = (6.66D-11)
  RCONST(1668) = (SUN*2.68E-4)
  RCONST(1669) = (1.4D-12*EXP(-1860/TEMP)*8.5)
  RCONST(1670) = (2.68D-11)
  RCONST(1671) = (SUN*1.99E-5)
  RCONST(1672) = (2.91D-13*EXP(1300/TEMP)*0.706)
  RCONST(1673) = (2.7D-12*EXP(360/TEMP))
! RCONST(1674) = constant rate coefficient
  RCONST(1675) = (6.00D-13*RO2)
  RCONST(1676) = (1.00D+06)
  RCONST(1677) = (2.91D-13*EXP(1300/TEMP)*0.706)
  RCONST(1678) = (2.7D-12*EXP(360/TEMP))
! RCONST(1679) = constant rate coefficient
  RCONST(1680) = (9.20D-14*RO2)
  RCONST(1681) = (1.4D-12*EXP(-1860/TEMP)*5.5)
  RCONST(1682) = (2.15D-11)
  RCONST(1683) = (SUN*1.99E-5)
  RCONST(1684) = (2.91D-13*EXP(1300/TEMP)*0.706)
  RCONST(1685) = (2.7D-12*EXP(360/TEMP))
! RCONST(1686) = constant rate coefficient
  RCONST(1687) = (9.20D-14*RO2*0.70)
  RCONST(1688) = (9.20D-14*RO2*0.30)
  RCONST(1689) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(1690) = (5.2D-13*EXP(980/TEMP)*0.56)
  RCONST(1691) = (7.5D-12*EXP(290/TEMP))
  RCONST(1692) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(1693) = (2.3E-12*1.74)
  RCONST(1694) = (1.00D-11*RO2)
  RCONST(1695) = (3.42D-11)
  RCONST(1696) = (SUN*5.79E-6+SUN*1.99E-5)
  RCONST(1697) = (SUN*5.79E-6)
  RCONST(1698) = (1.00D+06)
  RCONST(1699) = (2.91D-13*EXP(1300/TEMP)*0.625)
  RCONST(1700) = (2.7D-12*EXP(360/TEMP))
! RCONST(1701) = constant rate coefficient
  RCONST(1702) = (8.80D-13*0.2*RO2)
  RCONST(1703) = (8.80D-13*0.6*RO2)
  RCONST(1704) = (8.80D-13*0.2*RO2)
  RCONST(1705) = (6.99D-11)
  RCONST(1706) = (SUN*5.79E-6+SUN*5.50E-5)
  RCONST(1707) = (6.64D-11)
  RCONST(1708) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(1709) = (SUN*3.98E-6)
  RCONST(1710) = (SUN*5.79E-6)
  RCONST(1711) = (7.34D-12)
  RCONST(1712) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(1713) = (3.74D-12)
  RCONST(1714) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(1715) = (2.7D-12*EXP(360/TEMP))
! RCONST(1716) = constant rate coefficient
  RCONST(1717) = (8.00D-13*RO2)
  RCONST(1718) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(1719) = (2.7D-12*EXP(360/TEMP))
! RCONST(1720) = constant rate coefficient
  RCONST(1721) = (8.00D-13*RO2)
  RCONST(1722) = (1.53D-12)
  RCONST(1723) = (SUN*5.79E-6)
  RCONST(1724) = (7.54D-11)
  RCONST(1725) = (SUN*5.79E-6)
  RCONST(1726) = (1.00D+06)
  RCONST(1727) = (9.48D-11)
  RCONST(1728) = (SUN*5.79E-6)
  RCONST(1729) = (1.00D+06)
  RCONST(1730) = (9.13D-11)
  RCONST(1731) = (1.90D-12*EXP(190/TEMP))
  RCONST(1732) = (SUN*5.79E-6)
  RCONST(1733) = (1.00D+06)
  RCONST(1734) = (1.90D-12*EXP(190/TEMP))
  RCONST(1735) = (SUN*5.79E-6)
  RCONST(1736) = (1.00D+06)
  RCONST(1737) = (2.91D-13*EXP(1300/TEMP)*0.820)
  RCONST(1738) = (2.7D-12*EXP(360/TEMP))
! RCONST(1739) = constant rate coefficient
  RCONST(1740) = (9.20D-14*RO2)
  RCONST(1741) = (1.00D+06)
  RCONST(1742) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(1743) = (2.7D-12*EXP(360/TEMP))
! RCONST(1744) = constant rate coefficient
  RCONST(1745) = (8.00D-13*RO2)
  RCONST(1746) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(1747) = (2.7D-12*EXP(360/TEMP))
! RCONST(1748) = constant rate coefficient
  RCONST(1749) = (8.00D-13*RO2)
  RCONST(1750) = (1.53D-12)
  RCONST(1751) = (SUN*5.79E-6)
  RCONST(1752) = (7.54D-11)
  RCONST(1753) = (SUN*5.79E-6)
  RCONST(1754) = (1.00D+06)
  RCONST(1755) = (9.48D-11)
  RCONST(1756) = (SUN*5.79E-6)
  RCONST(1757) = (1.00D+06)
  RCONST(1758) = (9.13D-11)
  RCONST(1759) = (1.90D-12*EXP(190/TEMP))
  RCONST(1760) = (SUN*5.79E-6)
  RCONST(1761) = (1.00D+06)
  RCONST(1762) = (1.90D-12*EXP(190/TEMP))
  RCONST(1763) = (SUN*5.79E-6)
  RCONST(1764) = (1.00D+06)
  RCONST(1765) = (2.91D-13*EXP(1300/TEMP)*0.770)
  RCONST(1766) = (2.7D-12*EXP(360/TEMP))
! RCONST(1767) = constant rate coefficient
  RCONST(1768) = (8.00D-13*RO2)
  RCONST(1769) = (2.91D-13*EXP(1300/TEMP)*0.770)
  RCONST(1770) = (2.7D-12*EXP(360/TEMP))
! RCONST(1771) = constant rate coefficient
  RCONST(1772) = (8.00D-13*RO2)
  RCONST(1773) = (9.00D-13)
  RCONST(1774) = (SUN*5.79E-6)
  RCONST(1775) = (1.90D-12*EXP(190/TEMP))
  RCONST(1776) = (SUN*5.79E-6)
  RCONST(1777) = (1.00D+06)
  RCONST(1778) = (1.90D-12*EXP(190/TEMP))
  RCONST(1779) = (SUN*5.79E-6)
  RCONST(1780) = (1.00D+06)
  RCONST(1781) = (6.68D-11)
  RCONST(1782) = (SUN*5.79E-6)
  RCONST(1783) = (1.00D+06)
  RCONST(1784) = (1.23D-10)
  RCONST(1785) = (SUN*5.79E-6)
  RCONST(1786) = (1.00D+06)
  RCONST(1787) = (6.07D-11)
  RCONST(1788) = (9.18D-11)
  RCONST(1789) = (8.56D-11)
  RCONST(1790) = (SUN*5.79E-6)
  RCONST(1791) = (1.00D+06)
  RCONST(1792) = (1.42D-10)
  RCONST(1793) = (SUN*5.79E-6)
  RCONST(1794) = (1.00D+06)
  RCONST(1795) = (7.95D-11)
  RCONST(1796) = (1.11D-10)
  RCONST(1797) = (2.91D-13*EXP(1300/TEMP)*0.820)
  RCONST(1798) = (2.7D-12*EXP(360/TEMP))
! RCONST(1799) = constant rate coefficient
  RCONST(1800) = (8.00D-13*RO2)
  RCONST(1801) = (2.91D-13*EXP(1300/TEMP)*0.820)
  RCONST(1802) = (2.7D-12*EXP(360/TEMP))
! RCONST(1803) = constant rate coefficient
  RCONST(1804) = (8.00D-13*RO2)
  RCONST(1805) = (1.53D-12)
  RCONST(1806) = (SUN*5.79E-6)
  RCONST(1807) = (1.90D-12*EXP(190/TEMP))
  RCONST(1808) = (SUN*5.79E-6)
  RCONST(1809) = (1.00D+06)
  RCONST(1810) = (1.90D-12*EXP(190/TEMP))
  RCONST(1811) = (SUN*5.79E-6)
  RCONST(1812) = (1.00D+06)
  RCONST(1813) = (2.73D-12)
  RCONST(1814) = (6.19D-12)
  RCONST(1815) = (SUN*5.79E-6)
  RCONST(1816) = (1.12D-12)
  RCONST(1817) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(1818) = (1.00D+06)
  RCONST(1819) = (GCARR(2.66E-12,0.0E+00,200.0))
  RCONST(1820) = (GCARR(1.14E-12,0.0E+00,200.0))
  RCONST(1821) = (SUN*6.49E-6)
  RCONST(1822) = (GC_HACOHA(2.15E-12,0.0E+00,305.0))
  RCONST(1823) = (GC_HACOHB(2.15E-12,0.0E+00,305.0))
  RCONST(1824) = (SUN*2.88E-6)
  RCONST(1825) = (GCARR(4.00E-13,0.0E+00,0.0))
  RCONST(1826) = (1.88D-11)
  RCONST(1827) = (1.90D-12*EXP(190/TEMP))
  RCONST(1828) = (SUN*3.98E-6)
  RCONST(1829) = (SUN*5.79E-6)
  RCONST(1830) = (1.00D+06)
  RCONST(1831) = (5.86D-12)
  RCONST(1832) = (SUN*3.98E-6)
  RCONST(1833) = (1.72D-11)
  RCONST(1834) = (SUN*5.79E-6)
  RCONST(1835) = (1.00D+06)
  RCONST(1836) = (5.2D-13*EXP(980/TEMP)*0.56)
  RCONST(1837) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(1838) = (7.5D-12*EXP(290/TEMP))
  RCONST(1839) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(1840) = (2.3E-12*1.74)
  RCONST(1841) = (1.00D-11*RO2)
  RCONST(1842) = (1.90D-12*EXP(190/TEMP))
  RCONST(1843) = (SUN*5.79E-6)
  RCONST(1844) = (1.00D+06)
  RCONST(1845) = (1.90D-12*EXP(190/TEMP))
  RCONST(1846) = (SUN*5.79E-6)
  RCONST(1847) = (1.00D+06)
  RCONST(1848) = (1.15D-10)
  RCONST(1849) = (SUN*1.25E-4*2)
  RCONST(1850) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(1851) = (5.2D-13*EXP(980/TEMP)*0.56)
  RCONST(1852) = (7.5D-12*EXP(290/TEMP))
  RCONST(1853) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(1854) = (2.3E-12*1.74)
  RCONST(1855) = (1.00D-11*RO2)
  RCONST(1856) = (1.90D-12*EXP(190/TEMP))
  RCONST(1857) = (2.04D-13)
  RCONST(1858) = (SUN*2.68E-4)
  RCONST(1859) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(1860) = (5.2D-13*EXP(980/TEMP)*0.56)
  RCONST(1861) = (7.5D-12*EXP(290/TEMP))
  RCONST(1862) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(1863) = (2.3E-12*1.74)
  RCONST(1864) = (1.00D-11*RO2)
  RCONST(1865) = (5.2D-13*EXP(980/TEMP)*0.15)
  RCONST(1866) = (5.2D-13*EXP(980/TEMP)*0.41)
  RCONST(1867) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(1868) = (7.5D-12*EXP(290/TEMP))
  RCONST(1869) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(1870) = (2.3E-12*1.74)
  RCONST(1871) = (1.00D-11*RO2*0.3)
  RCONST(1872) = (1.00D-11*RO2*0.7)
  RCONST(1873) = (1.19D-11)
  RCONST(1874) = (SUN*5.79E-6)
  RCONST(1875) = (SUN*3.98E-6)
  RCONST(1876) = (1.00D+06)
  RCONST(1877) = (3.42D-11)
  RCONST(1878) = (SUN*5.79E-6)
  RCONST(1879) = (SUN*5.79E-6)
  RCONST(1880) = (SUN*1.99E-5)
  RCONST(1881) = (1.00D+06)
  RCONST(1882) = (5.2D-13*EXP(980/TEMP)*0.15)
  RCONST(1883) = (5.2D-13*EXP(980/TEMP)*0.41)
  RCONST(1884) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(1885) = (7.5D-12*EXP(290/TEMP))
  RCONST(1886) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(1887) = (2.3E-12*1.74)
  RCONST(1888) = (1.00D-11*RO2*0.3)
  RCONST(1889) = (1.00D-11*RO2*0.7)
  RCONST(1890) = (1.69D-11)
  RCONST(1891) = (SUN*5.79E-6)
  RCONST(1892) = (1.00D+06)
  RCONST(1893) = (1.34D-11)
  RCONST(1894) = (1.02D-13)
  RCONST(1895) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(1896) = (4.66D-11)
  RCONST(1897) = (SUN*5.79E-6)
  RCONST(1898) = (1.00D+06)
  RCONST(1899) = (2.55D-11)
  RCONST(1900) = (5.68D-12)
  RCONST(1901) = (1.90D-12*EXP(190/TEMP))
  RCONST(1902) = (SUN*5.79E-6)
  RCONST(1903) = (1.00D+06)
  RCONST(1904) = (1.90D-12*EXP(190/TEMP))
  RCONST(1905) = (SUN*5.79E-6)
  RCONST(1906) = (1.00D+06)
  RCONST(1907) = (7.04D-11)
  RCONST(1908) = (SUN*1.25E-4+SUN*2.68E-4)
  RCONST(1909) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(1910) = (5.2D-13*EXP(980/TEMP)*0.56)
  RCONST(1911) = (7.5D-12*EXP(290/TEMP))
  RCONST(1912) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(1913) = (2.3E-12*1.74)
  RCONST(1914) = (1.00D-11*RO2)
  RCONST(1915) = (1.52D-10)
  RCONST(1916) = (SUN*5.79E-6)
  RCONST(1917) = (1.00D+06)
  RCONST(1918) = (1.90D-12*EXP(190/TEMP))
  RCONST(1919) = (SUN*5.79E-6)
  RCONST(1920) = (1.00D+06)
  RCONST(1921) = (1.90D-12*EXP(190/TEMP))
  RCONST(1922) = (SUN*5.79E-6)
  RCONST(1923) = (1.00D+06)
  RCONST(1924) = (8.35D-11)
  RCONST(1925) = (SUN*1.25E-4+SUN*2.68E-4)
  RCONST(1926) = (5.2D-13*EXP(980/TEMP)*0.56)
  RCONST(1927) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(1928) = (7.5D-12*EXP(290/TEMP))
  RCONST(1929) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(1930) = (2.3E-12*1.74)
  RCONST(1931) = (1.0D-11*RO2)
  RCONST(1932) = (1.90D-12*EXP(190/TEMP))
  RCONST(1933) = (SUN*5.79E-6)
  RCONST(1934) = (1.00D+06)
  RCONST(1935) = (1.90D-12*EXP(190/TEMP))
  RCONST(1936) = (SUN*5.79E-6)
  RCONST(1937) = (1.00D+06)
  RCONST(1938) = (1.90D-12*EXP(190/TEMP))
  RCONST(1939) = (7.70D-11)
  RCONST(1940) = (SUN*1.25E-4*2)
  RCONST(1941) = (5.2D-13*EXP(980/TEMP)*0.56)
  RCONST(1942) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(1943) = (7.5D-12*EXP(290/TEMP))
  RCONST(1944) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(1945) = (2.3E-12*1.74)
  RCONST(1946) = (1.00D-11*RO2)
  RCONST(1947) = (9.58D-11)
  RCONST(1948) = (SUN*1.25E-4*2)
  RCONST(1949) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(1950) = (5.2D-13*EXP(980/TEMP)*0.56)
  RCONST(1951) = (7.5D-12*EXP(290/TEMP))
  RCONST(1952) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(1953) = (2.3E-12*1.74)
  RCONST(1954) = (1.00D-11*RO2)
  RCONST(1955) = (1.90D-12*EXP(190/TEMP))
  RCONST(1956) = (SUN*5.79E-6)
  RCONST(1957) = (1.00D+06)
  RCONST(1958) = (1.90D-12*EXP(190/TEMP))
  RCONST(1959) = (SUN*5.79E-6)
  RCONST(1960) = (1.00D+06)
  RCONST(1961) = (3.59D-12)
  RCONST(1962) = (SUN*5.79E-6)
  RCONST(1963) = (1.00D-14)
  RCONST(1964) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(1965) = (5.2D-13*EXP(980/TEMP)*0.56)
  RCONST(1966) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(1967) = (7.5D-12*EXP(290/TEMP))
  RCONST(1968) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(1969) = (2.3E-12*1.74)
  RCONST(1970) = (1.00D-11*RO2)
  RCONST(1971) = (9.83D-11)
  RCONST(1972) = (SUN*5.79E-6)
  RCONST(1973) = (9.51D-11)
  RCONST(1974) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(1975) = (1.00D+06)
  RCONST(1976) = (2.91D-13*EXP(1300/TEMP)*0.706)
  RCONST(1977) = (2.7D-12*EXP(360/TEMP))
! RCONST(1978) = constant rate coefficient
  RCONST(1979) = (2.00D-12*RO2)
  RCONST(1980) = (4.16D-12)
  RCONST(1981) = (SUN*5.79E-6+SUN*2.68E-4)
  RCONST(1982) = (5.66D-13)
  RCONST(1983) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(1984) = (6.58D-12)
  RCONST(1985) = (9.69D-12)
  RCONST(1986) = (SUN*5.79E-6)
  RCONST(1987) = (6.46D-12)
  RCONST(1988) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(1989) = (8.76D-13)
  RCONST(1990) = (4.34D-12)
  RCONST(1991) = (SUN*5.79E-6)
  RCONST(1992) = (7.36D-13)
  RCONST(1993) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(1994) = (5.2D-13*EXP(980/TEMP)*0.56)
  RCONST(1995) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(1996) = (7.5D-12*EXP(290/TEMP))
  RCONST(1997) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(1998) = (2.3E-12*1.74)
  RCONST(1999) = (1.00D-11*RO2)
  RCONST(2000) = (2.15D-11)
  RCONST(2001) = (SUN*3.98E-6+SUN*1.25E-4)
  RCONST(2002) = (2.75D-11)
  RCONST(2003) = (SUN*5.79E-6)
  RCONST(2004) = (2.42D-11)
  RCONST(2005) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(2006) = (5.2D-13*EXP(980/TEMP)*0.56)
  RCONST(2007) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(2008) = (7.5D-12*EXP(290/TEMP))
  RCONST(2009) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(2010) = (2.3E-12*1.74)
  RCONST(2011) = (1.00D-11*RO2)
  RCONST(2012) = (2.75D-11)
  RCONST(2013) = (SUN*5.79E-6)
  RCONST(2014) = (2.42D-11)
  RCONST(2015) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(2016) = (1.00D+06)
  RCONST(2017) = (5.2D-13*EXP(980/TEMP)*0.56)
  RCONST(2018) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(2019) = (7.5D-12*EXP(290/TEMP))
  RCONST(2020) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(2021) = (2.3E-12*1.74)
  RCONST(2022) = (1.00D-11*RO2)
  RCONST(2023) = (8.01D-11)
  RCONST(2024) = (SUN*5.79E-6)
  RCONST(2025) = (7.66D-11)
  RCONST(2026) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(2027) = (9.29D-11)
  RCONST(2028) = (SUN*5.79E-6)
  RCONST(2029) = (8.96D-11)
  RCONST(2030) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(2031) = (9.82D-11)
  RCONST(2032) = (4.15D-11)
  RCONST(2033) = (SUN*5.79E-6+SUN*1.25E-4)
  RCONST(2034) = (3.82D-11)
  RCONST(2035) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(2036) = (9.59D-12)
  RCONST(2037) = (SUN*5.79E-6+SUN*2.68E-4)
  RCONST(2038) = (1.00D+06)
  RCONST(2039) = (1.4D-12*EXP(-1860/TEMP)*4.0)
  RCONST(2040) = (1.27D-11)
  RCONST(2041) = (SUN*3.98E-6)
  RCONST(2042) = (3.60D-11)
  RCONST(2043) = (SUN*5.79E-6+SUN*1.25E-4)
  RCONST(2044) = (3.28D-11)
  RCONST(2045) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(2046) = (3.60D-11)
  RCONST(2047) = (SUN*5.79E-6+SUN*1.25E-4)
  RCONST(2048) = (3.28D-11)
  RCONST(2049) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(2050) = (3.06D-11)
  RCONST(2051) = (SUN*5.79E-6+SUN*1.25E-4)
  RCONST(2052) = (2.74D-11)
  RCONST(2053) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(2054) = (SUN*1.2e-5)
  RCONST(2055) = (SUN*1.2e-5)
  RCONST(2056) = (SUN*1.2e-5)
  RCONST(2057) = (SUN*1.2e-5)
  RCONST(2058) = (SUN*1.2e-5)
  RCONST(2059) = (SUN*1.2e-5)
  RCONST(2060) = (SUN*1.2e-5)
  RCONST(2061) = (SUN*1.2e-5)
  RCONST(2062) = (k_3rd(TEMP,2.5E+19_dp,9.7E-29,5.6,9.3E-12,1.5,0.6))
  RCONST(2063) = (k_3rd(TEMP,2.5E+19_dp,9.7E-29,5.6,9.3E-12,1.5,0.6))
  RCONST(2064) = (k_3rd(TEMP,2.5E+19_dp,9.7E-29,5.6,9.3E-12,1.5,0.6))
  RCONST(2065) = (k_3rd(TEMP,2.5E+19_dp,9.7E-29,5.6,9.3E-12,1.5,0.6))
  RCONST(2066) = (k_3rd(TEMP,2.5E+19_dp,9.7E-29,5.6,9.3E-12,1.5,0.6))
  RCONST(2067) = (k_3rd(TEMP,2.5E+19_dp,9.7E-29,5.6,9.3E-12,1.5,0.6))
  RCONST(2068) = (k_3rd(TEMP,2.5E+19_dp,9.7E-29,5.6,9.3E-12,1.5,0.6))
  RCONST(2069) = (k_3rd(TEMP,2.5E+19_dp,9.7E-29,5.6,9.3E-12,1.5,0.6))
  RCONST(2070) = (k_3rd(TEMP,2.5E+19_dp,9.7E-29,5.6,9.3E-12,1.5,0.6))
  RCONST(2071) = (k_3rd(TEMP,2.5E+19_dp,9.7E-29,5.6,9.3E-12,1.5,0.6))
  RCONST(2072) = (2.91D-13*EXP(1300/TEMP)*0.77*0.6)
  RCONST(2073) = (2.91D-13*EXP(1300/TEMP)*0.77*0.4)
  RCONST(2074) = (2.91D-13*EXP(1300/TEMP)*0.82*0.6)
  RCONST(2075) = (2.91D-13*EXP(1300/TEMP)*0.82*0.4)
  RCONST(2076) = (2.91D-13*EXP(1300/TEMP)*0.859*0.6)
  RCONST(2077) = (2.91D-13*EXP(1300/TEMP)*0.859*0.4)
  RCONST(2078) = (2.91D-13*EXP(1300/TEMP)*0.859*0.6)
  RCONST(2079) = (2.91D-13*EXP(1300/TEMP)*0.859*0.4)
  RCONST(2080) = (2.91D-13*EXP(1300/TEMP)*0.859*0.6)
  RCONST(2081) = (2.91D-13*EXP(1300/TEMP)*0.859*0.4)
      
END SUBROUTINE Update_RCONST

! End of Update_RCONST function
! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~


! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
! 
! Update_PHOTO - function to update photolytical rate constants
!   Arguments :
! 
! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

SUBROUTINE Update_PHOTO ( )


   USE aromatics_kpp_Global

      
END SUBROUTINE Update_PHOTO

! End of Update_PHOTO function
! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~



END MODULE aromatics_kpp_Rates

